BV SOLUTIONS AND VISCOSITY APPROXIMATIONS OF 
RATE-INDEPENDENT SYSTEMS 
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Abstract. In the nonconvex case solutions of rate-independent systems may develop jumps as 
a function of time. To model such jumps, we adopt the philosophy that rate independence should 
be considered as limit of systems with smaller and smaller viscosity. For the finite-dimensional 
case we study the vanishing-viscosity limit of doubly nonlinear equations given in terms of a 
differentiable energy functional and a dissipation potential which is a viscous regularization of 
a given rate-independent dissipation potential. 

The resulting definition of 'BV solutions' involves, in a nontrivial way, both the rate- 
independent and the viscous dissipation potential, which play a crucial role in the description 
of the associated jump trajectories. 

We shall prove a general convergence result for the time-continuous and for the time- 
discretized viscous approximations and establish various properties of the limiting BV solutions. 
In particular, we shall provide a careful description of the jumps and compare the new notion of 
solutions with the related concepts of energetic and local solutions to rate-independent systems. 
AMS Subject Classification; 49Q20, 58E99. 



1. Introduction 

Ratc-indcpcndcnt evolutions occur in several contexts. We refer the reader to [32| and the 
forthcoming monograph |39| for a survey of rate-independent modeling and analysis in a wide 
variety of applications, which may pertain to very different and far-apart branches of mechanics 
and physics. Rate-independent systems present very distinctive common features, because of 
their hysteretic character [54l [24] . Driven by external loadings on a time scale much slower than 
their internal scale, such systems respond to changes in the external actions invariantly for time- 
rescalings. Thus, they in fact show (almost) no intrinsic time-scale. This kind of behavior is 
encoded in the simplest, but still significant, example of rate- independent evolution, namely the 
doubly nonlinear differential inclusion 

a^olw'W) +D£t(wW) 3 inX* for a.a. t £ (0, T). (DNq) 

For the sake of simplicity, we will consider here the case when X is a finite dimensional linear 
space, £ : [Q,T] x X ^ M. an energy functional (D£ denoting the differential of £ with respect 
to the variable u e X), and "^q : X [0, -l-oo) is a convex, nondegcnerate, dissipation potential, 
hereafter supposed positively homogeneous of degree 1. Thus, ( |DNo[ ) is invariant for time-rescalings, 
rendering the system rate independence. 

Since the range K* of d'^o is a proper subset of X*, when £(t, •) is not strictly convex one 
cannot expect the existence of an absolutely continuous solution of ( |DNoD . Over the past decade, 
this fact has motivated the development of suitable notions of weak solutions to ( |DNoD . In the 
mainstream of [HI [35l [44], the present paper aims to contribute to this issue. Relying on the 
vanishing-viscosity approach, we shall propose the notion of BV solution to ( |DNo[ ) and thoroughly 
analyze it. 

To better motivate the use of vanishing viscosity and highlight the features of the concept of 
BV solution, in the next paragraphs we shall briefly recall the other main weak solvability notions 



Date: 9 October 2009. 

A.M. has been partially supported by DFG, Research Unit FOR 797 "MicroPlast" , Mi 459/5-1. R.R. and G.S. 
have been partially supported by a MIUR-PRIN'06 grant for the project "Variational methods in optimal mass 
transportation and in geometric measure theory" . 



2 



ALEXANDER MIELKE, RICCARDA ROSSI, AND GIUSEPPE SAVARE 



for dDNpl ). For the sake of simphcity, we shall focus on the particular case 

^o(w) ~ \\v\\, for some norm || ■ || on X. (1-1) 

Energetic and local solutions. The first attempt at a rigorous weak formulation of ( |DNoD goes 
back to [40j and the subsequent |421 141j , which advanced the notion of global energetic solution to 
the rate-independent system ( |DNo[ ). In the simplified case ()l.ip . this solution concept consists of 
the following relations, holding for all t G [0, T]: 

yzeX: Etiuit))<£t{z) + \\z-u{t)\\, (S) 

£t{u{t)) +Var{u;[0,t]) ^ EoHO)) + [ dt£s{u{s)) ds . (E) 

Jo 

The energy identity jEj balances at every time t G [0,T] the dissipated energy Var(M; [0,t]) (the 
latter symbol denotes the total variation of the solution u G BV{[0,T];X) on the interval [0,i]), 
with the stored energy Et{u{t)), the initial energy, and the work of the external forces. On the 
other hand, ([S| is a stability condition, for it asserts that the change from the current state u{t) to 
another state z brings about a gain of potential energy smaller than the dissipated energy. Since 
the competitors for u{t) range in the whole space X, ([S| is in fact a global stability condition. 

The global energetic formulation ©-jEl only involves the (assumcdly smooth) power of the 
external forces 9t£, and is otherwise derivative-free. Thus, it is well suited to jumping solutions. 
Furthermore, as shown in [571 [32], it is amenable to analysis in very general ambient spaces, even 
with no underlying linear structure. Because of its flexibility, this concept has been exploited 
in a variety of applicative contexts, like, for instance, shape memory alloys [HJ [371 Eji crack 
propagation [TH [Tl [T7] , elastoplasticity [Ml [301 [311 [201 [TOl El [2H1 , damage in brittle materials 
[38l [6l [521 [33] , delamination [23] , ferroelectricity [43] , and superconductivity (5^ . 

On the other hand, in the case of nonconvex energies condition ([S]) turns out to be a strong 
requirement, for it may lead the system to change instantaneously in a very drastic way, jumping 
into very far-apart energetic configurations (see, for instance, [30l Ex. 6.1], [211 Ex. 6.3], and [351 
Ex. 1]). On the discrete level, global stability is reflected in the global minimization scheme giving 
raise to approximate solutions by time-discretization. Indeed, for a fixed time-step r > 0, inducing 
a partition {0 = to < ti < . . . < t^-i < tj^ = T} of the interval [0,T], one constructs discrete 
solutions (U")^^i of ([S])-([E| by setting U° := uq and then solving recursively the variational 
incremental scheme 

ij;' G Argmin|||U-U;'-i|| + £f„(U)| for n = 1, . . . , iV. (IFq) 
vex ^ J 

However, a scheme based on local minimization would be preferable, both in view of numerical 
analysis and from a modeling perspective, see the discussions in |301 Sec. 6] and, in the realm of 
crack propagation, [16l [45l [26] . 

As pointed out in |16j . local minimization may be enforced by perturbing the variational 
scheme ( |IFo[ ) with a term, modulated by a viscosity parameter e, which penalizes the squared 
distance from the previous step U"~^ 

U;',GArgnmJ||U-U;'7'||+e^ + £*„ (U) [ for n = 1, . . . , iV , (IP^) 

vex I T J 

and depends on a second norm | ■ |, typically Hilbertian, on the space X. In a infinite-dimensional 
setting, one may think of X = L^(0), with a domain in W^, d> 1, and || • ||, | • | the and 
norms, respectively. Notice that, on the time-continuous level, ( |IFoD corresponds to the viscous 
doubly nonlinear equation 

^^^(^(O) +D£t(ue(t)) 3 in X* for a.a. t G (0,T), 

with ^^J^{v) = \\v\\ + ^-\vf ^ " 

(see [HI [5] for the existence of solutions G AC{[0,T]; X)). Then, the idea would be to consider 
the solutions to fDNoj) arising in the passage to the limit, in the discrete scheme PFgD , as e and 
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T tend to simultaneously, keeping t <C e. One can guess that, at least formally, this procedure 
should be equivalent to considering the limit of the solutions to ( |DNg[ ) as e | 0. 

Vanishing viscosity has by now become an established selection criterion for mechanically fea- 
sible weak solvability notions of rate-independent evolutions. We refer the reader to [25] for 
rate-independent problems with convex energies and discontinuous inputs, and, in more specific 
applied contexts, to [12] for elasto-plasticity with softening, to [H] for general material models with 
nonconvex elastic energies, the recent [13] for cam-clay non-associative plasticity, and [531 IHl 122] 
for crack propagation. Since the energy functionals involved in such applications are usually non- 
smooth and nonconvex, the passage to the limit mostly relies on lower semicontinuity arguments. 
Let us illustrate the latter in the prototypical case ( |DNg[ ). The key observation is that ( |DN^[ ) is 
equivalent (see the discussion in Section [2. 4p to the e-energy identity 

tt{uS))+ / (||<(s)llds+ J|<(s)p + -dist,(-D£,(Me(s)),if*)')ds 

= £o(u(0))+ / dt£.s{u,{s))ds 
Jo 

for all t G [0,T], where the term 

dist,(-D£t(u(i)),A'*) mill I -D£t(u(t)) -^1*, with K* = {z e X* : \\z\\, < l} , (1.3) 

measures the distance with respect to the dual norm | • j.^ of — D£t(u(i)) from the set K* . The term 
defined in (|1.3p is penalized in (|1.2p by the coefficient l/2e. Thus, passing to the limit in p.2p as 
e i 0, one finds 

dist*(-D£t(M(t)),X*) = for a.a.t e (0,r) . 

Hence, 

-D£t(u(0) G i^*, i.e. II -D£t(u(i))||, < 1 for a.a. t G (0, T) , (1.4) 

which is a local version of the global stability Furthermore, p.2p yields, via lower-semicontinuity, 
the energy inequality 

£t(u(t)) + Var(u; [0,t]) < £o(u(0)) / dt£.s{u{s)) ds for alH G [0, T] . (1.5) 

Jo 

Conditions (|1.4p - (|1.5p give raise to the notion of local solution of the rate-independent sys- 
tem ( jPN^ . 

While the local stability (|1.4p is more physically realistic than jS]), its combination with the 
energy inequality (|1.5p turns out to provide an unsatisfactory description of the solution at jumps 
(see the discussion in [35l See. 5.2] and Remark 12.81 later on). In order to capture the jump 
dynamics, the energetic behavior of the system in a jump regime has to be revealed. From this 
perspective, it seems to be crucial to recover from p.2p . as e ], 0, an energy identity, rather than 
an energy inequality. Thus, the passage to the limit has to somehow keep track of the limit of the 
term 




which in fact encodes the contribution of the viscous dissipation ||u^p, completely missing in (|1.5p . 
BV solutions. Moving from these considerations, it is natural to introduce the vanishing viscosity 
contact potential (which is related to the bipotential discussed in [7] , see Section [3]) induced by 
^e, i.e. the quantity 




= ||w|| + \v\ dist^{w,K*) tor V e X , w e X* . 



Then, the £-energy identity (|1.2I) yields the inequality 

E.tiue{t))+ f p«(s),-D£,(M,(s)))ds<£o(u(0))+ / 9t£.(?/e(s))ds , (1.7) 
Jo Jo 
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see Section [3.11 Passing to the limit in p.7p . in Theorem 14.101 we shall prove that, up to a 
subsequence, the solutions (ite) of the viscous equation ( |DNg[ ) converge, as e J, 0, to a curve 
u e BV([0, r];X) satisfying the local stability (|1.4p and the following energy inequality 

£t(u(t)) + yarp.£(u;[0,t])<£oKO))+ / 9t£,,(w(s)) ds . (1.8) 

Jo 

Without going into details (see Definition 13.41 later on), we may point out that (|1.8p features a 
notion of (pscudo)-total variation (denoted by Varp ^) induced by the vanishing viscosity contact 
potential p (|1.6|1 and the energy £. The main novelty is that a BV-curve obeying the local stability 
condition ()1.4p always satisfies the opposite inequality in (jl.Sp . thus yielding the energy balance 

£t(ii(t)) + yarp.£(u;[t,t])==£o(u(0))+ / dtt s{u(s)) ds . (Ep,£) 

Jo 

In fact, Varp_£ provides a finer description of the dissipation Ap.£ of u, along a jump between 
two values u_ and at time t: it involves not only the quantity — it_|| related to the 
dissipation potential (jl.ip . but also the viscous contribution induced by the vanishing viscosity 
contact potential p through the formula 

Ap,£(t; u^,u+) := inf | H p(79(r), -D£t(79(r))) dr : 

^•^'■0 (1.9) 

^ e AC([ro,ri];X), d{ro) = m_, i9(ri) = «+}• 

By a rescaling technique, it is possible to show that, in a jump point, the system may switch to a 
viscous behavior, which is in fact reminiscent of the viscous approximation ( |DNg[ ). In particular, 
when the jump point is of viscous type, the infimum in (|1.9p is attained and the states U- and 
M+ are connected by some transition curve •& : [ro, ri] — > X, fulfilling the viscous doubly nonlinear 
equation 

a«'o(t5'(r)) +i9'(r) +D£t(z?(r)) 9 in X* for a.a. r G (ro, n) 

(in the case the norm | • | is Euclidean and we use its differential to identify X with X*). The combi- 
nation of (|1.4p and (|1.8p yields the notion of BV solution to the rate-independent system {X, £, p). 
This concept was first introduced in [35| . in the case the ambient space X is a finite-dimensional 
manifold X, and both the rate-independent and the viscous approximating dissipations depend 
on one single Finsler distance on X. In this paper, while keeping to a Banach framework, we 
shall considerably broaden the class of rate- independent and viscous dissipation functionals, cf. 
Remark 12.41 Moreover, the notion of BV solution shall be presented here in a more compact form 
than in [35], amenable to a finer analysis and, hopefully, to further generalizations. 

Let us now briefly comment on our main results. First of all, we are going to show in Theo- 
rems HUl 1121 and 14. 71 that the concept of BV rate- independent evolution completely encompasses 
the solution behavior in both a purely rate- independent, non-jumping regime, and in jump regimes, 
where the competition between dry-friction and viscous effects is highlighted. Indeed, from (|1.4p 
and (jl.Sp it is possible to deduce suitable energy balances at jumps (cf. conditions ( |JbvP in 
Theorem ESI). 

Then, in Theorem 14. 101 we shall prove that, along a subsequence, the viscous approximations 
arising from dmQ converge as e: J, to a BV solution. Next, our second main result. Theo- 
rem |3?TT1 states that, up to a subsequence, also the discrete solutions Ur.e constructed via the 
e-discretization scheme dlP^P converge to a BV solution u S BV([0,T];X) of ( |DNoP as e | and 
T I simultaneously, provided that the respective convergence rates are such that 

£ 

lim — = +00 . 

Finally, in Section Owe shall develop a different approach to BV solutions, via the rescaling 
technique advanced in [18| and refined in |35[ I44j . The main idea is to suitably reparametrize the 
approximate viscous curves (u^) in order to capture, in the vanishing viscosity limit, the viscous 
transition paths at jumps points. This leads to performing an asymptotic analysis as e J, of 
the graphs of the functions u^, in the extended phase space [0,r] x X. For every e > the 
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graph of can be parametrized by a couple of functions (t^, Ug), t^ being the (strictly increasing) 
rescaling function and := o t^ the rescaled solution. In Theorem 15.61 we assert that, up to a 
subsequence, the functions (t^, Ue) converge as e | to a parametrized rate-independent solution. 
By the latter terminology we mean a curve (t, u) : [0, S] [0, T] x X fulfilling 

t : [0, S] [0, T] is nondecreasing. 



t'(s) + ||u'(s)|| > for a.a.s G (0,S), 



(1.10a) 



^ ||-D£,(,)(u(s))|| <1, 1 foraasGfOS) fl 10b) 

||u'(.)||>0 ^ ||-D£,(,)(u(.))||>l / fo'^a.a.sG(0,5), (1.10b) 

and the energy identity 

A£(t(s),u(s))-9,£(t(.),u(s))t'(s) ^^^^^^ 

= -||u'(s)|| - |u'(s)|dist,(-D£t(,)(u(s)),/s:*) for a.a. s e (0, S) , 

As already pointed out in [T51 [3S] , like the notion of B V solution, relations p.lOp as well comprise 
both the purely rate-independent evolution as well as the viscous transient regime at jumps. The 
latter regime in fact corresponds to the case — D£t(u) ^ K* : the system does not obey the local 
stability constraint (|1.4p any longer, and switches to viscous behavior, see also Remark 15.71 later 
on. 

As a matter of fact. Theorem 15.81 shows that parametrized rate-independent solutions may be 
viewed as the "continuous counterpart" to BV evolutions. With a suitable transformation, it 
is possible to associate with every parametrized rate-independent solution a BV one, and con- 
versely. One advantage of the parametrized notion is that it avoids the technicalities related to 
BV functions. Hence, it is for instance more easily amenable to a stability analysis (cf. [35l Rmk. 
6]). Furthermore, in [44j a highly refined vanishing viscosity analysis has been developed, with 
this reparametrization technique, in the infinite-dimensional (L^, L^)-framework, where ( |DNgP is 
replaced by a general quasilinear evolutionary PDE. 

Generalizations and future developments. So far we have focused on dissipation functionals 
of the type (|l.ip and ^'e(f) = ||w|| + as in dPN^D for expository reasons only, in order to 

highlight the main variational argument leading to the notion of BV solution. Indeed, the analysis 
developed in this paper is targeted to a general 

positively 1-homogeneous, convex dissipation Vfo : X — ^ [0, -l-oo), 

(cf. (|2.1[) '). and considers a fairly wide class of approximate viscous dissipation functionals ^'g, 
defined by conditions (j'I'.ip - (j'I'.3[) in Section 12.31 Furthermore, at the price of just technical 
complications, our results could be extended to the case of a Finsler-like family of dissipation 
functionals '^q{u, •), depending on the state variable u S X, and satisfying uniform bounds and 
Mosco-eontinuity with respect to u, see [351 Sect. 2] and [371 Sect. 6, 8]. 

The extension to infinite- dimensional ambient spaces and nonsmooth energies is crucial for ap- 
plication of the concept of BV solution to the PDE systems modelling rate-independent evolutions 
in continuum mechanics. A first step in this direction is to generalize the known existence results 
for doubly nonlinear equations, driven by a viscous dissipation, to nonconvex and nonsmooth en- 
ergy functionals in infinite dimensions. As shown in [481 147j . in the nonsmooth and nonconvex 
case one can replace the energy differential D£t with a suitable notion of subdifferential dEf Ac- 
cordingly, instead of continuity of D£(, one asks for closedness of the multivalued subdifferential 
dEt in the sense of graphs. These ideas shall be further advanced in the forthcoming work [55] . 
Therein, exploiting techniques from nonsmooth analysis, we shall also tackle energies which do not 
depend smoothly on time (this is relevant for rate-independent applications, see e.g. [H] and ^E^). 

On the other hand, the requirement that the ambient space is finite-dimensional could be 
replaced by suitable compactness (of the sublevels of the energy) and reflexivity assumptions 
on the ambient space X. The latter topological requirement in fact ensures that X has the 
so-called Radon-Nikodym property, i.e. that absolutely continuous curves with values in X are 
almost everywhere differentiable. The vanishing viscosity analysis in spaces which do not enjoy 
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this property requires a subtler approach, involving metric arguments (see e.g. [471 Sect. 7]), or 
ad-hoc stronger estimates [33]. See also [53] for some preliminary approaches to BV solutions for 
PDE problems. 

Plan of the paper. Section [2] is devoted to an extended presentation of energetic and local solu- 
tions to rate-independent systems. In particular, after fixing the setup of the paper in Scction [2.11 
in Sec. l2.2l we recall the definition of global energetic solution, show its differential characterization 
and the related variational time-incremental scheme. We develop the vanishing-viscosity approach 
in Sees. 12.31 and 12.41 thus arriving at the notion of local solution (see Section r2.5|l . which also 
admits a differential characterization. 

In Section [3] we introduce the concept of vanishing viscosity contact potential and thoroughly 
analyze its properties, as well as the induced (pseudo)-total variation. With these ingredients, in 
Sec.|4]we present the notion of BV solution. We show that BV rate-independent evolutions admit, 
too, a differential characterization, and, in Sec. 14.21 that they provide a careful description of the 
energetic behavior of the system. Then, in Section [4. 3( we state our main results on BV solutions. 

While SectionOis focused on the alternative notion of parametrized rate-independent solutions, 
the last Sec. [S] contains some technical results which lie at the core of our theory. 

2. Global energetic versus local solutions, and their viscous regularizations 

In this section, we will briefly recall the notion of energetic solutions and show that their viscous 
regularizations give raise to local solutions. 

2.1. Rate-independent setting: dissipation and energy functionals. We let 

{X, II • ||x) be a finite-dimensional normed vector space, 
endowed with a gauge function 'I'o, namely a 

non-degenerate, positively 1-homogeneous, convex dissipation ^'o : ^ ^ [0, -l-oo), (2.1) 
i.e. \['o satisfies 4*0 («) > if 7^ 0, and 

folwi +V2) < 'i'oivi) + *o(i'2), *o(Aw) = A^'o(w) for every A > 0, v,vi,V2 e X. 
In particular, there exists a constant 77 > such that 

jy^'^ll^'ILY < ^'o(i') < '7l|v||x for every v £ X. 
Since is 1-homogeneous, its subdifferential d'^o : X =^ X* can be characterized by 

d-^oiv) |ui e X : {w, z) < ^0(2) for every z e X, {w,v) = 4'o(v)} C X*; (2.2) 

d'^Q takes its values in the convex set K* C X*, given by 

/\* = 9*0(0) := {w e X* : (^,2) < l'o(z) Vz e X} D a^'o(w) for every w e X, (2.3) 

which enjoys some useful (and well-known, see e.g. [46j ) properties. For the reader's convenience 
we list them here: 

Kl. K* is the proper domain of the Legendre transform of 5*0, since 

X , 1 if w e K*, 

%{w)^lK'iw)^{ . (2.4) 

I +00 otlierwise. 

K2. 'So is the support function of K*, since 

4*0 (w) = sup {w,v) for every w g X, (2.5) 

weK* 

and K* is the polar set of the unit ball if := |w G X : 'I'o(u) < l} associated with ^Pq. 
K3. K* is the unit ball of the support function ^fo, of K: 

A'* = |w G X* : «'o,(w) < 1), with ^(^.(u;) = sup(w,w) = sup Vtt- (2-6) 
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K4. In the even case (i.e., when ^o(w) = ^'o(^w) for aU v G X), we have that is an equivalent 

norm for X, 5'o» is its dual norm. K and K* are their respective unit balls. 
Further, we consider a smooth energy functional 

£ G Ci([0,T] X X), 

which we suppose bounded from below and with energy-bounded time derivative 

3C> V(t,u) G [0,T] X X : £t(w)>-C, |9t£t(w)| < C (l + £t(w)+) , (2.7) 

where (•)+ denotes the positive part. The rate-independent system associated with the energy 
functional £ and the dissipation potential can be formally described by the rate- independent 
doubly nonlinear differential inclusion 

a*o(w'(i)) +D£t(u(t)) 9 inX* for a.a. t G (0, T). (DNq) 

As already mentioned in the Introduction, for nonconvex energies solutions to ( |DNo[ ) may 
exhibit discontinuities in time. The first weak solvability notion for ( |DNo[ ) is the concept of (global) 
energetic solution to the rate-independent system ( |DNo[ ) (see [3111301111] and the survey [32]), 
which we recall in the next section. 

2.2. Energetic solutions and variational incremental scheme. 

Definition 2.1 (Energetic solution). A curve u G BV([0, T\]X) is an energetic solution of the rate 
independent system (X, £,5'o) if for all t G [0,r] the global stability (S) and the energy balance 
(E) holds: 

VzgX: £t(M(t)) < £t(z) + ^o(2-w(i)), (S) 
£t(M(t))+Var^„(u;[0,t]) = £o(u(0))+ / dt£- s{u{s)) As . (E) 

BV functions. Hereafter, we shall consider functions of bounded variation pointwise defined in 
every point t G [0,T], such that the pointwise total variation with respect to 4*0 (any equivalent 
norm of X can be chosen) Var 152(11; [0,T]) is finite, where 

Var*o(ii; [a,b]) := sup | ^ *o(w(im) - M(tm-i)) : a = io < < ' ' ' < ^M-i < tn = &|- 

m—l 

Notice that a function u in BV{[0,T]; X) admits left and right limits at every t G [0,T] : 
u{t^) := limu(s), u(i+) := limu(s), with the convention u(0_) := u(0), u{T^) :~ u{T), (2.8) 

s\t sit 

and its pointwise jump set is the at most countable set defined by 

J„ :== {t G [0,r] : u{t^) ^ u{t) or u{t) ^ u{t+)] Z) ess-J„ {t G [0,T] : w(i_) ^ u{t+)]. (2.9) 

We denote by u' the distributional derivative of w, and recall that u' is a Radon vector measure 
with finite total variation \u'\. It is well known that u' can be decomposed into the sum of the 
three mutually singular measures 

u' = u'c^ + u'q + u'j, u'c^ = iiJ^^ ^ := u'^ + u'q . (2.10) 

Here, u'^ is the absolutely continuous part with respect to the Lebesgue measure , whose 
Lebesgue density u is the usual pointwise (and ^-^-a.e. defined) derivative, Mj is a discrete measure 
concentrated on ess-J„ C Ju, and u'q is the so-called Cantor part, still satisfying it^di}) = for 
every t G [0,r]. Therefore u'^^ = u'^ + u'q is the diffuse part of the measure, which does not 
charge J„. In the following, it will be useful to use a nonnegative and diffuse reference measure ^ 
on (0,T) such that and u'q are absolutely continuous w.r.t. ^: just to fix our ideas, we set 

^i:=^' + \u'q\. (2.11) 
With a slight abuse of notation, for every (a, 6) C (0,T) we denote by /^^d4'o(uco) the integral 

C.^/*.Wd^V/'*.(||) d|„fc|, (2.2, 
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Since ^I^o is 1-homogcneous, the above integral is independent of /i, provided u'^^ is absolutely 
continuous w.r.t. /i. 

Towards a differential characterization of energetic solutions. Let us first of all point out 
that ([S| is stronger than the local stability condition 

- D£t(w(i)) e K* for every t G [0, T] \ J„, (Sioc) 

which can be formally deduced from ( |DNo[ ) and ()2.3|) . Indeed, the global stability ([S| yields for 
every z — u{t) + hv G X and h > Q 

(-D£t(u(t)), hv) + o{\h\) < £t(u(t)) - £t(u(t) + hv) < h^o{v) 

and therefore, dividing by h and passing to the limit as /i | 0, one gets 

{-B£t{u{t)),v) < -^oiv) for every z e X, 

so that dSiocD holds. We obtain more insight into (|E| by representing the 'I'o variation Var^p^ (u; [a, b]) 
in terms of the distributional derivative u' of u. In fact, recalling (|2.1ip and (|2.12p . we have 

Varvi,„(u; [a,6]) := / d*o«o) + JmP*o("; [«' ^])' 

J a 

where the jump contribution Jmp,p^^(u; [a, b]) can be described, in terms of the quantities 

A*o(i'o,wi) ^o(«i - Wo), A^,,{v_,v,v+) ■=-9oiv~v_) + '^oiv+-v), (2.13) 

by 

Jmp^„(M;[a,6]) A*„(u(a),u(a+)) + A^.„(w(6_),w(6))+ ^ A^,„{u{t_),u{t),u{t+)). (2.14) 

tGJun(a,6) 

Also notice that, as usual in rate- independent evolutionary problems, u is pointwise everywhere 
defined and the jump term Jmp^^(u; [•,•]) takes into account the value of u at every time t € J^^. 
Therefore, if u is not continuous at t, this part may yield a strictly bigger contribution than 
the total mass of the distributional jump measure Uj (which gives rise to the so-called essential 
variation) . 

The following result provides an equivalent characterization of energetic solutions: besides the 
global stability condition ([5]), it involves a BV formulation of the differential inclusion ( |DNoD (cf. 
the subdifferential formulation of [41j ) and a jump condition at any jump point of u. 

Proposition 2.2. A curve u S BV([0,r];X) satisfying the global stability condition jS]) is an 
energetic solution of the rate-independent system {X, £, \l/o) if and only if it satisfies the differential 
inclusion 

dfi 

and the jump conditions 
£t{u{t))-£t{u{t^)) = -A^,{u{t_),u{t)), £t(u(i+))-£t(u(t)) = -A>t„(u(t),u(i+)), 

£.t{u{t+)) - ^t{u{t^ j) = -AMt-),Ht+))- ^ 
for every t G Ju (recall convention (j2.8p in the case t = 0,T). 

We shall simply sketch the proof, referring to the arguments for the forthcoming Proposition 12.71 
for all details. 

Proof. By the additivity property of the total variation Yar^g{u; [■,■]), jEj yields for every < 
to<ti<T 

Var*„(«;[io,ii]) + £ti(«(ii)) = £t„(w(to))+ dtEt{u{t)) At . (E') 

J to 

Arguing as in the proof of Proposition 12. 71 later on, one can see that the global stability ^ and 
(|E^ yield the differential inclusion (DNq.bv) and conditions dJcncrl )- 



o(^(i)) +D£t(u(t)) 9 for ^i-a.e.t&[Q,T], ^Ji := + \u'c\, (DNo,bv) 



Conversely, repeating the arguments of Proposition [^771 one can verify that (DNq^bv I and dJcnciD 



imply dE]). □ 



BV SOLUTIONS AND VISCOSITY APPROXIMATIONS OF RATE-INDEPENDENT SYSTEMS 



9 



Incremental minimization scheme. Existence of energetic solutions can be proved by solving 
a minimization scheme, which is also interesting as construction of an effective approximation of 
the solutions. 

For a given time-step r > we consider a uniform partition (for simplicity) = to < < ■ ■ ■ < 
tAr_i < T < In, tn := ?^T, of the time interval [0,r], and an initial value w uq. In order to 
find good approximations of U" ~ u(tn) wc solve the incremental minimization scheme 

find U^, • • ■ , such that U;' G Argmin |*o(U - U;'"^) + £t„ (U)|. (IPo) 

vex ^ J 

Setting 

Vr{t):=V'^ if f e (t„_i,t„], (2.15) 

it is possible to find a suitable vanishing sequence of step sizes r^. | (see, e.g., [JU [32] for all 
calculations), such that 

3 lim Ur,(t) =: w(t) for every t e [0, T], 

k — *+oo 

and u is an energetic solution of ( |DNo[ ). 

2.3. Viscous approximations of rate-independent systems. In the present paper we want to 
study a different approach to approximate and solve ( |DNo[ ) : the main idea is to replace the linearly 
growing dissipation potential with a suitable convex and superlinear "viscous" regularization 
^'e : X ^ [0, +00) of ^0, depending on a "small" parameter e > and "converging" to in a 
suitable sense as e | 0. Solving the doubly nonlinear differential inclusion (we use the notation u 
for the time derivative when u is absolutely continuous) 

a*e(we(i)) + D£t(we(t)) 3 inX* for a.a. t e (0, T), (DN^) 

one can consider the sequence (u^) as a good approximation of the solution u of ( |DNo[ ) as e | 0. 

There is also a natural discrete counterpart to ( |DNg[ ), which regularizes the incremental mini- 
mization problem ( |IPo[ ) . We simply substitute ^'o by in ( |IPo[ ) , recalling that now the time-step 
T should explicitly appear, since ^'e is not 1-homogeneous any longer. The viscous incremental 
problem is therefore 

find U^.,--- ,Uf^ such that U^J ^ G Argmin (t*^ f^^^^ — + £*„ (U)}. (IP^) 

vex I \ T / J 

Setting as in 

\JrAt)-^Vle if i e (<„-l,t„], 

one can study the limit of the discrete solutions when r j and e i 0, under some restriction on 
the behavior of the quotient e/r (see Theorem 14 . 1 1 1 later on). 

The choice of the viscosity approximation ^E'g. Here we consider the particular case when 
the potential can be obtained starting from a given 

*(«) 

convex function ^ : X ^ [0, +00) such that *(0) = 0, hm 77-17^ = +00, ('I'.l) 

ll"IUT+oo \\v\\x 

by the canonical rescaling 

^'j(w) e~i5i(eti) for every v G X, e > 0, (5'. 2) 

and is linked to 'I'o by the relation 

vI/q(v) = lim\I'£(w) = lime^^^{ev) for every v e X. (\E'.3) 

Remark 2.3. Notice that, by convexity of and the fact that ^'(0) = 0, the map e 1— > e~^^'(ew) 
is nondecreasing for all v £ X. Hence, 

«'o(w) < ^'^(u) for all w G a:, for all £ > 0. (2.16) 

Furthermore, by the coercivity condition (j^.ip . 

d'^s{v) := 94' (eu) is a surjective map. 
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Here are some examples, showing that (|^.2[) stih provides a great flexibihty and covers several 
interesting cases. 

Example 2.4. 

v&o-viscosity: The simplest example, still absolutely non trivial [35], is to consider 



d^e{v) = (l + e*o(«))5*o(«). 



(2.17) 



A similar regularization can be obtained by choosing a real convex and superlinear function 
Fv : [0, +oo) [0, +oo), with Fv{0) = i^y(O) = 0, and setting 

^{v) := *o(w) + Fv{-fo{v)) = F(*o(w)), with F{r) := r + Fv{r). (2.18) 

Quadratic or p- viscosity induced by a norm || • ||: The most interesting case involves an ar- 
bitrary norm || • || on X and considers for p > 1 

^iv)^^o{v) + -\\v¥, *s(«) = *o(t') + ^^r, d^,{v) = d-^o{v) Jp{v), (2.19) 
P P 

where Jp is the p-duality map associated with || ■ ||. In particular, if || • || is a Hilbertian norm 
and p ~ 2, then J2 is the Riesz isomorphism and we can choose J2iv) ~ v by identifying X 
with X*. Hence, JmQ reads 

d^,{u,{t)) + eu,{t) +J:)£.t{u,{t)) 3 in X* for a.a. t G (0,T), 

and the incremental problem ( |IPe[ ) looks for J7"^ which recursively minimizes 

u ^ Mu - K]-') + ^\\u- + ^tAu). 

This is the typical situation which motivates our investigation. 
Additive viscosity: More generally, we can choose a convex "viscous" potential ^'y : X 
[0, -l-cxo) satisfying 

lime~i*y(£u) = 0, lim A"i*y(Av) = +00 for aU v eX, (2.20) 

eiO Af+cx) 

and set 

-^{v) ■.= ^o{v) +'^v{v), ^,{v) ■.^■9oiv)+e-^^v{ev), ^^'^(w) = 9*0 + (9*y (ew). (2.21) 

2.4. Viscous energy identity. Since "if has a superlinear growth, the results of j51[S] ensure that 
for every e > and initial datum uq ^ X there exists at least one solution E AC{[0,T]; X) to 
equation fDN^]), fulfilhng the Cauchy condition u^{0) = uq. 

In order to capture its asymptotic behavior as e | 0, we split equation ( |DNgD in a simple 
system of two conditions, involving an auxiliary variable Wg : [0, T] —>■ X* and a scalar function 
p, : [0,T] 

d-^,{u,{t)) 3 w, for a.a. t G (0, T) , (2.22a) 
mt{u,{t)) ^ ^w,it), dttt{ue{t)) = -Pe{t) for alH G [0,r]. (2.22b) 

Denoting by ^'*, ^* the conjugate functions of 'I' and ^f^, we have 

= **(0) < **(0 < +00, **(0=e"^**(0 for every e € (2.23) 

Due to (f2l6l) . there holds 

*e(0<*o(0 for all G a:, e > 0. (2.24) 

The classical characterization of the subdifferential of ^'e yields that the first condition (|2.22ap is 
equivalent to 

■^e{ue{t)) + K{we{t)) = {we{t),ue{t)) for a.a. t G (0, T) . (2.25) 
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On the other hand, the chain rule for the functional £ shows that along the absolutely contin- 
uous curve Me 

^£t(Me(i)) = {BEt{u,it)),u,{t))+dtet{ue{t)) = -{w,{t),u,{t))-p,{t) for a.a. t e (0,T). (2.26) 
Thus, if Welt) ~ — D£f(ue(t)), equation (|2.22ap is equivalent to the energy identity 

(u,(r)) + VI/* {w,{r))+p,{r)) dr + S-tAMh)) = E-toMo)), (2.27) 

for every < to < ^i < T"- 

Remark 2.5 (The role of v]/*). In the general, additive-viscosity case (see (|2.2ip ). when 4'(f) = 
vE'o(w) + vl>y(w) the inf-sup convolution formula yields 

*:(0= inf (l^.(a) + ivI/^(6)) mmM/t.(e-z). 

In particular, when '^viO ■= sI'^^P some norm | • | of X, one finds 

ze zsA * 

where | • |* is the dual norm of | • |. Thus, for all ^ G X* the functional vl'*(^) is the squared distance 
of ^ from K*, with respect to | • |*. This shows that, in the viscous regularized equation ( |DNgD , 
the (local) stability condition w{t) = — D£((it(i)) e K* has been replaced by the contribution of 
the penalizing term 

^/ min |~D£t(u,(t))-z|2dt 

in the energy identity (|2.27p . 

2.5. Pointwise limit of viscous approximations and local solutions. Using (|2.7p . it is not 

difficult to show that the viscous solutions of ( |DNgP satisfy the a priori bound 

^ («'e(ue(t)) + «'*(we(0)) dt < C, with (t) = -DE,t{ue{t)) for all i G [0,r]. (2.28) 

Therefore, Helly's compactness theorem shows that, up to the extraction of a suitable subsequence, 
the sequence (u^) pointwise converges to a BV curve u. From the convergence Wg{t) — > w{t) = 
-D£t(M(i)) as e i and the fact that for aU t G [0, T] 



STm if n)(t) G K* 

limMe-'^*{we{t)) > %{w{t))^rMt)) = \ ' (2.29) 

eio I +00 otherwise, 

we infer that the limit curve u satisfies the (local) stability condition ( |Sioc[ ). On the other hand, 
passing to the limit in (|2.27p one gets the energy inequality 

£t,(ii(ti))+Var*„(u;[to,ii]) < £io("(io))+ / dt£t{u{t))dt for < to < ^i < T. (E;„,q) 

Jto 

The above discussion motivates the concept of local solution (see also j35[ Sec. 5.2] and the refer- 
ences therein). 

Definition 2.6 (Local solutions). A curve u G JiV{[0,T]; X) is called a local solution of the rate 
independent system {X, £, 5'o) if it satisfies the local stability condition 

- mt{u{t)) G K* for every t G [0, T] \ J„ (Sioc) 



and the energy dissipation inequality (E| 



mcq 



Local solutions admit the following differential characterization. 
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d^o{^{t)) +DE.t{u{t))3 forf,-a.e.te[0,T], /i if^ + (DNo,bv) 



Proposition 2.7 (Differential characterization of local solutions). A curve u G BV([0, T]]X) is a 
local solution of the rate independent system (X, £, 4'o) if and only if it satisfies the BV differential 
inclusion 

d^ 

and the jump inequalities 
£t{u{t)) - £t{u{t-)) < -A^„{u{t^),u{t)), Et{u{t+)) - EtHt)) < -A^„{u{t),u{t+)), 
et{u{t+)) - < -A^,{u{t^),u{t+)), 

at each jump time i G Ju • 



(Jlocal) 



Proof. Notice that at every point t € (0,r) where du'^^{t) / dfi = 0, the differential inclusion 
(DNq.bvI reduces to the local stability condition dSiod ). In the general case, (DNq^bv) follows by 
differentiation of (£(^^^,^1. Indeed, the latter procedure provides the following inequality between 
the distributional derivative ^£t(u(t)) of the map t £t{u{t)) and the 'fp-total variation measure 
*o«o) := *o(d<JdM)/^ for ^ := u'^ + 



-£t(M(t)) + *o«J - dtlt{u{t))S£^ < . 



(2.30) 



Applying the chain rule formula for the composition of the functional £ and the BV curve u 
(see [2] and [3l Thm. 3.96]) and taking into account the fact that u'^^ and Uj are mutually singular, 
we obtain from (|2.30p that 



(2.31) 



we finally conclude (DNq bvI- Localizing (£1^^^,^) around a jump point t we get 



Combining (|2.3ip with the local stability condition fSj^, in view of the characterization (12.21) of 
9*0 and of ([! 
the inequalities dJiocail ) 

Conversely, let us suppose that a BV curve u satisfies (DNq^bv I and dJiocaiD - The local stability 
condition is an immediate consequence of (DNq.bvK which yields — D£t(u(i)) G K* for 
a.e. t G [0, T] and therefore, by continuity, at every point of [0, T] \ J„. 

In order to get (Einoq'i again apply the chain rule for the composition £ and u, obtaining 



£ti (u(ii)) + 



-B£t{u{t)) 



d/i 



d/i(t) - Jmp(£; [to,ii]) 



(2.32) 



£to(«(^o))+ / dtEtHt)) dt, 



where 



and 



Jmp(£;[<o,ti]) = S+(to) + -B_(ii) 



E 

tGJun(to,ti) 



E-{t) 



E_{t) := Et{u{t)) - £t(u(i-)), E+{t) := £*(«(<+)) - £.t{u{t)). 



By (DNo,BV I we have 

(-D£,(«(t)), df,{t) = 1^'^ *o(^(t)) dM(t) = d*0«o) : 

whereas dJiocail ) yields for every t G 

E^{t) < -A^„{u{t-),u{t)), E+{t) < ^A^,{u{t),u{t+)), 



so that — Jmp(£; [toi^i]) > Jmp*o(w; [toj^i]) and therefore (£(^^^^1 follows from (|2.32p . 



(2.33) 

(2.34) 
□ 
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Remark 2.8. Unlike the ease of energetic solutions (cf. Proposition 12. 2p . a precise description of 
the behavior of local solutions at jumps in missing here. In fact, the jump inequalities dJiocail ) are 
not sufficient to get an energy balance and do not completely capture the jump dynamics, see the 
discussion of [35l Sec. 5.2]. 

In order to get more precise insight into the jump properties and to understand the correct 
energy balance along them, we have to introduce a finer description of the dissipation. It is 
related to an extra contribution to the jump part of Var^(,(it; [•,•]), which can be better described 
by using the vanishing viscosity contact potential induced by the coupling 4", 5"*. We describe this 
notion in the next section. 

3. VANISHING VISCOSITY CONTACT POTENTIALS AND FiNSLER DISSIPATION COSTS 

3.1. Heuristics for the concept of vanishing viscosity contact potential. Suppose for the 
moment being that, in a given time interval [ro,ri], the energy £t(-) = £(•) does not change w.r.t. 
time. If "i? e AC{[ro,ri]; X) is a solution of ( |DNgD connecting uq = i?(ro) to ui = i?(ri), then the 
energy release between the initial and the final state is, by the energy identity (|2.27p . 

£K)-£(ui)=/ + *;H)dt, 

Jro ^ ' (3.1) 

with i>(^) = iSl(t) and = -D£(i9(t)) for a.a. t G (0, T). 

If one looks for a lower bound of the right-hand side in the above energy identity which is inde- 
pendent o/e > 0, it is natural to recur to the functional p : X x X* [0, +oo) defined by 

p(?;, w) := inf (^'^(w) + ^%{w)) = inf (e~^^'(ew) + £"^^'*(w)) for u e X, w G X*. 
We obtain 

£(mo)-£(wi)> / p(u,w)di with = and = -D£(t9(t)). (3.2) 

J ro 

Since p(-,-) is positively 1-homogeneous with respect to its first variable, the right-hand side 
expression in (j3.2p is in fact independent of (monotone) time rescalings. On the other hand, the 
vanishing viscosity contact potential p(-, •) has the remarkable properties 

\){v,w) > {w,v), \){v,w) >^o{v) for every w e X, w g X*. (3.3) 

Therefore, ii-d G AC{[rQ,ri];X) is another arbitrary curve connecting mq to ui, the chain rule (|2.26p 
for £ yields 

£{uo)-£{ui)= f \w{t),v{t))dt< f \'fMt))+'^*e{wm dt 

Jro Jro 

(where v denotes the time derivative of and w = — D£(i?)), whence 

£(uo) - £(wi) < / p{v{t),w{t))dt. (3.4) 

Jro 

It follows that, in a time regime in which the energy functional £ does not change with respect to 
time, for every £ > any viscous solution of ( |DNg[ ) (and, therefore, any suitable limit of viscous 
solutions) should attain the minimum dissipation, measured in terms of the vanishing viscosity 
contact potential p. Moreover, this dissipation always provides an upper bound for the energy 
release, reached exactly along viscous curves and their limits. 

Remark 3.1. In some of the cases discussed in Example 12.41 the vanishing viscosity contact po- 
tential p admits a more explicit representation. 

(1) We first consider the "if o -viscosity case (|2.18p . where ^(u) := F{'i'o{v)), F : [0, -f-cx)) 
[0,-|-oo) being a real convex superlinear function with F{0) = 0,F'{0) = 1. We introduce the 
1-homogeneous support function ^Pq, of the set 

K := {v eX : *o(w) < l}, 4'o*(w;) := sup{w,v). 

veK 
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It is not difficult to show that 'i'*{w) = F*(*o*(w')) and that for all {v, w) e X x X* 

p(«,«;) = *oMmax(l,M'o,M) = , , ![ ^ (3.5) 

[*o(w)*o*(w) ifw^JsT*. 

(2) In the additive viscosity case of (|2.2ip one has for all {v, w) G X x X* 

p{v,w) ~ '^q{v) + pv{v,w), where pv{v,w) = inf { e^^'i'viev) + e^^ inf '^y{w — z)]. (3.6) 

In particular, when "^viv) ~ Fv{\\v\\) for some norm || • j| of X and a real convex and superlinear 
function Fy : [0, +oo) [0, +oo) with (0) = F{,(0) = 0, we have for aU {v, w) e X x X* 

p(u, w) = ^'o(w) + pv(«j w), with pv(z;,w) = ||u|| min Hw — z||,. (3.7) 



Notice that in (|3.5p and (j3.7p the form of the vanishing viscosity contact potential p does not 
depend on the choice of F and Fy, respectively, but only on the chosen viscosity norm. 

By the 1-homogcneity of p(-, w) and these variational properties, it is then natural to introduce 
the following Finsler dissipation. 

Definition 3.2 (Finsler dissipation). For a fixed t G [0, T], the Finsler cost induced by p and (the 
differential of) £ at the time t is given by 

Ap.£(<;uo,wi) :=inf I / ' p(^9(r), -D£t(^(r))) dr : 

^ -^-o (3.8) 

z9e AC([ro,ri];X), z?(ro) = uq, ^n) = u,^ 
for every wo,wi G X. We also consider the induced "triple" cost 

Remark 3.3. Since p(v,w) > 5'o(w) by p.3p . a simple time rescaling argument shows that the 
infimum in (|3.8p is always attained by a Lipschitz curve G AC([ro, ri]; X) with constant p- 
speed, in particular such that 

p(iil(r),-D£t(i?(r)) = 1 for a.a. r G (ro,ri) . 

By the heuristical discussion developed throughout (|3.ip - (|3.4p . the cost Ap.£ is the natural 
candidate to substituting the potential ^E'o and the related cost A^^, of (|2.13p in the jump con- 
tributions (|2.14p and in the jump conditions dJcncrP - Notice that the second relation of (|3.3p 
implies 

Ap,£(t;wo,tti) > A*o(mo,'"i) for every Mo, Ml G X. (3.9) 
The notion of jump variation arising from such replacements is precisely stated as follows. 

Definition 3.4 (The total variation induced by Ap_£). Let u G BV([0, T];X) a given curve, let 
be the diffuse part of its distributional derivative u' , and let be its pointwise jump set (j2.9p . 
For every subinterval [a, b] C [0, T] the Jump variation of u induced by (p, £) on [a, b] is 

Jmpp £(u; [a, 6]) :=Ap_£ (a; u(a), 'u(a+)) + Ap,£ (6; u(6_), -u(6))-|- 

+ \Mt;u{t.),u{t),uit+)), (3-10) 

teJun(a,6) 

and the (pseudo-)total variation induced by (p, £) is 

Varp,E{u;[a,b]):^ f d1'o«o) + JmPp,£ (^i [a, &])• (3.11) 

J a 

Remark 3.5 (The (pseudo-)total variation Varp^£). Let us mention that Varp_£ enjoys some of 
the properties of the usual total variation functionals, but it is not lower semicontinuous w.r.t. 
pointwise convergence. In fact, it is not difficult to see that its lower semicontinuous envelope is 
simply Var^Q. Furthermore, Varp^£ is not induced by any distance on X. Indeed, we have used 
slanted fonts in the notation Var to stress this fact. In order to recover a more standard total 
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variation in a metric setting, one has to work in the extended space !^ :— [0, T\y. X and add 
the local stability constraint — D£t G K* on the "continuous" part of the trajectories. We shall 
discuss this point of view in Section [6l 

In view of inequality p.9p between the Finsler dissipation Ap^£ and Aip^, the notion of total 
variation associated with Ap,£ provides an upper bound for Var^^, namely 

Vm e BV([0,T];X), [a,6]c[0,r]: Varp^s (m; [a, 6]) > Var*,, (u; [a, 6]). (3.12) 

3.2. Vanishing viscosity contact potentials. While postponing the definition of BV solutions 
related to Varp^g to the next section, let us add a few remarks about the vanishing viscosity contact 
potential p 

p(w,w) inf (^-^fw) + = inf (e-i^few) +e"i**H) for u £ X, w <^ X* . (3.13) 

which partly matches the definition introduced by [7] . We first list a set of intrinsic properties of 
p, which we shall prove at the end of this section. 

Theorem 3.6 (Intrinsic properties of p). The continuous functional p : XxX* [0, +oo) defined 
by p.l3p satisfies the following properties: 

(11) For every v £ X,w £ X* the maps p(u, ■) and p(-, w) have convex sublevels. 

(12) p(w, w) > {w, v) for every v & X,w G X* . 

(13) For every w S X* the map v i— > p(y,w) is 1-homogeneous and thus convex in X, with 
p(v,w) >Oifv^O. 

(14) For every v £ X,w £ X* the map X i— *■ p{v, Xw) is nondecreasing in [0, +oo). 

(15) If for some vq £ X and iD, w G X* we have p{vo, w) < p{vo,w), then the inequality p{v, w) < 
p{v,w) holds for every v £ X, and there exists vi £ X such that p{vi,w) < {w,vi). 

Remark 3.7 (A dual family of convex sets). Property (15) has a dual geometric counterpart: let 
us first observe that for every w £ X* the map v ^ p{v, w) is a gauge function and therefore it is 
the support function of the convex set 

:= |z G X* : {z,v) < p{v,w) for every v £ x|, i.e. p{v,w) ~ sup|(z,w) : z £ ii'^j. 

Assertion (15) then says that for every couple w,w £ X 

we always have w £ K^, or w £ A'*, and, moreover, w £ AT*, <^ p(-,w) < Pi'jw). (3.14) 

Suppose in fact that w ^ K*^: this means that an element £ X exists such that {w,vq) > 
p{vo,w); by (12) we get p{vq,uj) > p{vo,w), and therefore by (15) p{v,w) > p{v,w) > {w,v) for 
every v £ X, so that w £ AT*,. The second statement of p.l4p is an immediate consequence of the 
second part of (15). 

Property (12) suggests that the set where equality holds in plays a crucial role: 
Definition 3.8 (Contact set). The contact set Sp C A x X* is defined as 

Ep := |(w,w) G A X A* : p{v,w) = (3.15) 

Here are some other useful consequences of (I1~I5) 

Lemma 3.9. If p : X x X* ^ [0, +oo) satisfies (I1~I5), then 

(16) for every v £ X,w £ X* we have 

p{v, 0) + Iks (^) ^ P(«' ^) ^ P("' 0)- (3-16) 

(17) The contact set can be characterized by 

(u,w)Gl]p w£dpi-,w){v) <^ v£dlKi{w). (3-17) 

More generally, if w £ dp{-,w){v) then {v,w) £ Sp, w £ K^, and p(v,w) = p{v,w). In 
particular, if w £ dK^ then w £ K*^ . 
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Proof. The chain of incquaUties in (|3.16p is an immediate consequence of (14) and of (|3.14p . (|3.17|1 
is a direct consequence of the fact that v >—> p{v,w) is a gauge function and Ik^ is its Legendre 
transform. 

In order to check the last statement, given w S X, u; G X* let us take w G dp{-,w){v) so that 
w G and p{v,w) — {iu,v). Combining (12) with ()3.14p we get p{v,w) = p{v,w), so that 
(w, w) G Sp. □ 

Remark 3.10. Properties (11,12,15) suggest a strong analogy between p and the notion of bipotential 
introduced by [7j: according to [7], a bipotential is a functional b : X x X* (— oo, +oo] which 
is convex and lower semicontinuous in each argument, satisfies (12), and whose contact set fulfils 
a condition similar to (|3.17p 

{v,w)eT,b ^ w e db{-,w){v) ^ V e db{v,-){w). 

In our situation, (|3.17|1 is a direct consequence of the homogeneity of p, but the convexity condition 
with respect to w looks too restrictive, as shown by this simple example. Consider the case 
X = X* = with *(v) := |[w||i + ^'y(u), ||u||i := \vi\ + jval, and 
11 13 

By ([3?6l) we have p{v,w) = \\v\\i+pv{v,w) with py(w,w) = infe>o 7 (*y(ew)+**(w) ) and find 



**(^)-^(KI-i)^ + ^(KI-i)f . 

Considering the special case v = (vi, 0), w ~ (0, W2), we obtain 

2/3 



py((z;i,0),(0,W2)) - v/372|«i|((k'2| -1)+)^ 

The map 1^2 > p((ui,0), {0,W2)) is therefore not convex. 

Let us now consider some properties of p and its contact set Ep involving explicitly the func- 
tional 4*. Since the vanishing viscosity contact potential p is defined through the minimum proce- 
dure (|3.13p . the contact set is strictly related to the set of optimal e > attaining the minimum 

in (ixra . 

Definition 3.11 (Lagrange multipliers). For every [v, w) € X x X* we introduce the multivalued 
function A (with possibly empty values) 

A{v,w) := je > : p(u,w) = «'e(w) + *e(w)| C [0,+oo). (3.18) 

Notice that for every {v,w) £ X x X* the function e > e~^'^{ev) + e^^'^*{w) is convex on 
(0, -l-oo). Since ^E" has superlinear growth at infinity, it goes to -l-oo as e t +00 if w ^ 0, so that 

the set A{v, w) is always a bounded closed interval if w ^ 0. (3.19) 
Tiieorem 3.12 (Properties of p, and Sp). 

(PI) The vanishing viscosity contact potential p satisfies p{v,0) — ^'o(t'), Kq = K* , and in 
particular 

p{v,w) > {w,v), ^o{v) + Ik'-{w) >p{v,w) >'^o{v) >0 for every v G X, w G X* , (3.20) 

p{v,w) ^-^oiv) ^ weK*. (3.21) 
(P2) For every w G X*, the convex sets are the sublevels of 

= jz G X* : **(z) < **(w)}, (3.22) 

and p admits the dual representation 

p(w,u;) = sup|(z,u) : z G X*, **(z) < **(w)|. (3.23) 

In particular, '^*{wi) < 'I'*(w2) for some wi, W2 G X* if and only i/p(w, wi) < p(v,W2) for 
every v € X. 
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(P3) The multivalued function A defined in (j3.18p is upper semicontinuous, i.e. 

if {vn,Wn) —>■ {v,w) € X X X* flTirf £„ S A(i;„ , w„) — > e, then eeA(w,w). (3.24) 

(P4) The contact set Sp (|3.15p can be characterized by 

m; G 5*o(w) C A'* or, if K* , 3e > : w e d^^ev), (3.25) 

and the last inclusion holds exactly for e G A{v,w). Equivalently, 

(v,w) e Sp <S4> e (95'e(u) for every e G K(v,w). 

In particular, in the case of additive viscosity, with'^{v) ~ '^q{v) + '^v{v) and ^'y satisfying 
(|2.20p . we simply have 

(w, w) e Ep 3 A > : we (9*o(w) + <9*y(Aw). (3.26) 
Proo/s of Theorems\SIMand\EIE 

Ad (PI). Inequalities p.20[) are immediate consequences of the definition of p. The equal- 
ity ^'o(f) = p{v,w) is equivalent to the existence of a sequence > such that (recall that 

lim £^^^'(£fct;) = l'o(w), lim £^^**(w) = 0. 

— 'oc k — 'oc 

Since the first inequality prevents Sk from diverging to +00 (being 5" superlinear), from the second 
limit we get ^'*(w) = 0, i.e. 

(u;,z)<*(z) VzGA:. 

Replacing z with ez, multiplying the previous inequality by e~^, and passing to the limit as £ J, 0, 
in view of (j^.3|) we conclude 

(«;,z)<*q(z) VzGX, so that w £ K* . 

The converse implication in (|3.2ip is immediate. 

Ad (P2). Since the sublevcls of ^P* are closed and convex, a duality argument shows that (|3.22p 
is equivalent to (|3.23p . In order to prove the latter formula, let us observe that, if 5'*(z) < 5'*(u'), 
then (2, v) < p{v, w), because the Fenchel inequality yields 

{z,v) = £"^(z,£w) < £~^^'(£?;) +£~^5'(z) = 5'e(u) + ^'*(z) < *e(u) + ^'e(u') for every £ > 0. 

We show that there exists z e X* such that ^*{z) < '^'*{w) and p{v,w) = {z,v). Due to (|3.2ip . if 
w G K* , then p{v, w) = 5'o(^^) and the thesis follows from (|2.5p Hence, let us suppose that w ^ K* 
and V ^ 0] then we can choose £0 G h.{v, w), Sq > 0, such that 

p{v,w) = £0 ^*(£ow) + £o"^^'*(w) < £"^^'(£v) + £"^^'*(it)) for every £ > 0. (3.27) 

Choosing z^ G d'^{ev) we have 

4'(£u) — 5'(£ov) < (^e, (£ — £o)w) for every £ > 

so that, in view of inequality p.27p . 

(e^^ - £(7^) (^^'(£ow) + 1'*(w)^ +£"^(ze, (£ - £o)u) > for every £ > 0. 

Dividing by £ — £o and passing to the limit first as £ J, £o and then as £ | £o, we thus find 
z± G 95'(£ou) (accumulation points of the sequences {z^ : e > Eq) and {z^ : e < £o), respectively), 
such that 

{z_,v) <piv,w)^eo'[^ieav) + ^*iw)^ < {z+,v). (3.28) 
On the other hand, the Fenchel identity of convex analysis yields 

£q ^**(z) = (z, v) - £(7^*(£ow) for every z G d^{sov) (3.29) 
so that the map z i— > 'i*{z) is affine on d'i'{eQv) and a comparison between p.28p and p.29p yields 

^*{z^) < **(w) < **(z+). 
Using formula (|3.29p we can thus find 9 G [0, 1] and zg := (1 — 6)z^ + Oz+ G (95'(£ow) such that 
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The last statement of (P2) follows easily. One implication is immediate. On the other hand, if 
> 'i'*{w2); then by the Hahn-Banach separation theorem we can find v G X and 5 > 
such that 

{wi,v) >S + {z,v) for every z S X* such that 5'*(z) < 'i'*{w2), 

and. therefore, by (|3.23p we conclude \){v,wi) > {wi,v) > 5 + p{v,W2). 
Ad (11,2,3,4,5) These properties directly follow from (P2). 

Ad (P3) and continuity of p. Notice that p is upper semicontinuous, being defined as the 
infimum of a family of continuous functions. Take now converging sequences (t^n), (w^n); (^n) as 
in (|3.24p : we have that 



liminf (e,/*(£„w„) (w„)) > 'f,{v) + ^:{w) 

n — >oo 

Since 



e-^*(£u) +£"^**(w) if£>0, 
*o(w) if£ = 0. 



p{v,w) > hmmf p(u„,u;„) > limmf (e„ *(£„u„) + £„ 

«^oo «^oc. (3 30) 

> + %iw) > lp{v,w), 

we obtain e G A(v, w). 

Inequality (|3.30p shows that p is also lower semicontinuous, since, if u ^ 0, any sequence £„ G 
A(w„,?ii„) admits a converging subsequence, in view of (|3.19p . 

Ad (P4). Concerning the characterization p.25p of Sp, it is easy to check that, if {v,w) satisfies 
p.25[) . then by the Fenchel identity and formula (|2.2[) we have, when w G K*, 

p(u, w) > {w, v) = *o(w) = p{v, w), 

and, when w ^ A'*, 

p{v,w) > (w,v) = e^^(w,ev) = £"^*(£v) +e^^^*{w) > p(w,w) 

so that {v,w) G Sp and e G A{v,w). Conversely, if p{v,w) — {w,v) and w G K*, then by (|3.20p 
'i>o{v) = {w, v) and therefore w G d'i/o{v). If w ^ K*, then, choosing e G A(u, w), we have 

^(£w) + = ep{v, w) = {w, ev), so that w G 9^(£w). 

In the particular case of (|2.2ip , (|3.26p follows now from (|3.25p by the sum rule of the subdifferentials 
and the 0- homogeneity of 95'o- Q 



4. BV SOLUTIONS AND ENERGY-DRIVEN DISSIPATION 

4.1. BV solutions. We can now give our precise definition of BV solution of the rate-independent 
system (X, £,p), driven by the vanishing viscosity contact potential p (|3.13p and the energy £. 
From a formal point of view, the definition simply replaces the global stability condition jS]) by 
the local one dSiod ), and the ^'o-total variation in the energy balance ([E| by the "Finsler" total 
variation ()3.1ip . induced by p and £. 

Definition 4.1 (BV solutions, variational characterization). A curve u G BV([0,T];X) is a BV 
solution of the rate independent system (A", £,p) the local stability dSiocP and the (p, E)-energy 
balance hold: 

~mt{u{t))eK* fora.a.fG [0,T]\ J„ (Sioc) 

yarp.£(u;[0,t]) + £t(w(t)) = £o(u(0))+ / 9t£,,(u(s)) ds for alH G [0, T]. (Ep.g) 

Jo 

We shall see in the next Section l473l that any pointwise limit, as £ | 0, of the solutions (u^) of the 
viscous equation ( |DNeP or, as r, £ J, 0, of the discrete solutions (Ur.e) of the viscous incremental 
problems PP^P , is a BV solutions induced by the vanishing viscosity contact potential p. Let us 
first get more insight into Definition 14. II 
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Properties of BV solutions. As in the case of energetic solutions, it is not difficult to see that 
the energy balance ( |Ep g ) holds on any subinterval [^0,^1] C [0,r]; moreover, if the local stability 
condition dSiod ) holds, to check (Ep_£ ) it is sufficient to prove the corresponding inequality. 



Proposition 4.2. // u G BV([0,T];X) satisfies (Ep g), then for every subinterval [to,ti] there 
holds 

yarp,£(u;[<o,ii]) + £ti(u(ti)) = £to(M(io))+ / dtEsHs)) ds. (E; 

Jt„ 

Moreover, if u satisfies ( |SiocD , then ( Ep^£[ ) is equivalent to the energy inequality 

yarp,£(u; [0,T]) + £.t{u{T)) < EoHO)) + f dt£.s{u{s)) ds. (Ep,£;ineq) 



Proof. ( Ep £ I easily follows from the additivity property 

VO < to < ii < i2 < T : yarp,£(u; [to,ti]) + Varp,£(u; [^1,^2]) = Var^x{u; [^0,^2]). (4.1) 

In order to prove the second inequality we argue as in |351 Prop. 4], taking dSiod ) into account. □ 

Notice that, by (|3.12p . any BV solution is also a local solution according to Definition l2.61 i.e. it 
satisfies the local stability condition and energy inequality (Ei,^^. In fact, one has a more accurate 



description of the jump conditions, as the following Theorem shows (cf. with Propositions 12.21 
andEIZl). 



Theorem 4.3 (Differential characterization of BV solutions). A curve u G BV([0, T]\X) is a BV 
solution of the rate-independent system (X, £,p) if and only if it satisfies the doubly nonlinear 
differential inclusion in the BV sense 

fdu' 



-(t)) +D£t(u(t)) 9 for li-a.e. t e[Q,T], := ^ + lu'd. 



and the following jump conditions at each point t G J,i of the jump set (|2.9I 

£,t{u{t)) - £t(w(t-)) = -Ap,£(i;u(t_),u(t)), 
^t{u{t+)) - £-t{u{t)) = -Ap,£(t;ii(t),u(t+)), 
£t(?/(i+))-£t(u(t_)) = -Ap,£(t;u(t_),u(t+)). 



(DNo.Bv) 



(Jbv) 



Proof. We have already seen (see Lemma [2.7p that local solutions satisfy (DNq^bvI- The jump 
conditions ( |JbvD can be obtained by localizing (E^ ^ ) around any jump time t G 3u- 

Conversely, to prove (Ep £.i,ioq I (as seen in the proof of Lemma l2.7[ dSiod ) ensues from (DNq.bv 0; 
we argue as in the second part of the proof of Lemma 12.71 still applying (|2.32p and (|2.33p . but 
replacing inequalities (|2.34p with the following identities. 



E-it) 



(t.;u{t^),u{t)), E+{t) ^ -Ap,E{t;u{t),u{t+)) for alH G 



which are due to ( |JbvI )- Hence, --Jmp(£; [0, T]) = Jmpp £(u; [0, T]). Then, (Ep £.ineq I follows from 

m^- ' ' ' □ 

The next section is devoted to a refined description of the behavior of a BV solution along the 
jumps. 

4.2. Jumps and optimal transitions. Let us first introduce the notion of optimal transition. 

Definition 4.4. Let t G [0, T], G X with -D£t(it_), -D£f(u+) G A'*, and -00 < ro < 

ri < +00. An absolutely continuous curve iD : [ro,ri] —^ X connecting w_ = i9(ro) and it+ ~ t?(?'i) 
is an optimal (p, Et) -transition between u_ and it+ if 



d{r)^0 fora.a. rG(ro,ri); «'o*(-D£t(t9(r))) > 1 VrG[ro,ri], 



£((w_) - £t(it+) = Ap,£(t;M_,u+) 



p(i9(r),-D£t(i?(r)))dr. 



(0.1) 
(0.2) 
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Wc also say that an optimal transition ?9 is of 

sliding type if - D£t(i9(r)) G K* for every r g [ro, ri], (OsUding) 

OTscoits type if - D£t(?9(r)) ^ A'* for every r G (rg , ri ) , (Oviscous) 

energetic type if £t(u+) - £t(u_) = -'J'o(t^+ - "-)■ (Ocncr) 

We denote by 0(t;u_,M+) the (possibly empty) collection of such optimal transitions, with nor- 
malized domain [0, 1] and constant Finsler velocity 

p{d{r),-Btt{'d{r)))=Et{u-)~S.t{u+) for a.a. r G (0, 1) . (4.2) 

Remark 4.5. Notice that the notion of optimal transition is invariant by absolutely continuous 
(monotone) time rescalings with absolutely continuous inverse; moreover, any optimal transition 'd 
has finite length, it admits a reparametrization with constant Finsler velocity p(i?(-)7 ~D£t('i9(-))), 
and is a minimizer of (|3.8p . so that it is not restrictive to assume d G 0(t, u_,u+). 



Theorem 4.6. A local solution u G BV([0,T];X) is a BV solution according to Definition 
if and only if at every jump time t G the initial and final values u{t^) and u{t^) can be 
connected by an optimal transition curve i?* G Q{t]u{t-),u(t^)), and there exists r G [0,1] such 
that u{t) = t?*(r). Any optimal transition curve i} satisfies the contact condition 

(z9(r), -D£t(i9(r))) G Sp for a.a. r G (0, 1). (4.3) 

Proof. Taking into account Theorem l4.3[ the proof of the first part of the statement is immediate. 
To prove ()4.3p . let t be a jump point of u and let us first suppose that u{t-) = u(t) ^ u(i+)- 
By Remark 13.31 can find a Lipschitz curve i?oi G AC([ro, ri]; X) with normalized speed 
p(i?, — D£t('!9)) = 1, connecting w(t-) to M(i+), so that the jump condition ( |JbvD yields 

\-B£tmr))J{r))dr = £t(u(t_)) - £t{u{t+)) = /" ' p(79(r), -D£t(79(r))) dr. 

ro J ro 

This shows that § is an optimal transition curve and satisfies 

^p(i?, -D£t(i?(r))) dr - (-D£t(z?(r)), d{r))) dr = 0. 



Since the integrand is always nonnegative, it follows that (|4.3p holds. 

In the general case, when u is not left or right continuous at t, we join two (suitably rescaled) 
optimal transition curves i?oi G 9(t; u(t_), u(i)) and di2 G Q{t\u{t),u{tj^)). □ 

The next result provides a careful description of (p, £f)-optimal transitions. 

Theorem 4.7. Let t G [0,T], 7i_, u+ G X, and i? : [0, 1] ^ X he an optimal transition curve in 
Q{t;u^,u+). Then, 

(1) § is a constant- speed minimal geodesic for the (possibly asymmetric) Finsler cost Ap_£ (t; u_, u+), 
and for every < po < pi < 1 it satisfies 

e-tMpo)) - ^tmpi)) = A,.i,{t;d{po),^{pi)) 

f \ (4-4) 

= {Pi - Po)Ap^£(t;it_,M+) = [pi - p^)[Et{u-) - et{u+)j; 

In particular, the map p i-^ £t(i?(/?)) is affine. 



(2) An optimal transition is of sliding type (Ogiiding) if md only if it satisfies 

d-^oi^ir)) +Y)Et{'&{r)) 3^ /or a.a. r G (0, 1), (4.5) 

*o*(-D£i(i?(?'))) = 1 for every r G [0, 1]. (4.6) 

(3) An optimal transition -d is of viscous type dOyiscousT ) */ o-nd only if there holds for every selection 
(0, 1) 9 r e{r) in A(i$l(r), -D£t(i?(r)) 

5^'(£(r)tj(r)) +D£t(i?(r)) 3 /or a.a. r G (0, 1). (4.7) 
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Equivalently, there exists an absolutely continuous, surjective time resettling r : (po,Pi) 
(0, 1), with — oo < po < pi < oo and r{s) > for a.e. s S (poiPi); such that the rescaled 
transition 6{s) :~ 'd{r{s)) satisfies the viscous differential inclusion 

9*(0(s)) +D£t(6l(s)) 3 for a.a.s e {po,pi) , with lim6'(s)=M_, \im 9{s)^u+. (4.8) 

slpo sTPl 

(4) Any optimal transition 'd can be decomposed in a canonical way into an (at most) countable 
collection of optimal sliding and viscous transitions. In other words, there exists (uniquely 
determined) disjoint open intervals {S.j)j!^a and {Vk)k£v of (0,1), with cr, u C N, such that 
(0, 1) C ( \Jjea 'S'j) U ( Ufce^, Vk) and 

is of sliding type, '&\y^ 'is of viscous type. 

(5) An optimal transition is of energetic type dOcnciD */ cind only if d is of sliding type and it 
is a '^Q-minimal geodesic, i.e. 

*o(i?(''i) - -Oiro)) = (ri - ro)1'o(ui - u^) for every < ro < ri < 1. (4.9) 

//^'o has strictly convex sublevels, then § is linear and r {''^^{r)-, ^ti'&i'f))) is a linear segment 

contained in the graph of Et. 

If is Gateaux- differentiable at X \ {0} then 

-D£t(i?(r)) = D*o(^t+ - U-) for every r e [0, 1]. 

In particular, the map r t-^ — D£t(i9(r)) is constant. 

Remark 4.8. It follows from the characterization in (2) of Thcorcni l4.7l (cf. with (|4.5p ~ (|4.6p ) that 
sliding optimal transitions are independent of the form of the vanishing viscosity contact potential 
p, and thus on the particular viscosity potential 

Instead, as one may expect, occurs in the doubly nonlinear equation (|4.7p (equivalently, 
in (|4.8p ). which in fact describes the viscous transient regime. Hence, different choices of the 
viscous dissipation ^ shall give raise to a different behavior in the viscous jumping regime, see 
also the example in [5TJ Sec. 2.2]. The latter paper sets forth a different characterization of 
rate-independent evolution, still oriented towards local stability, but derived from a global-in-time 
variational principle and not a vanishing viscosity approach. 

Proof. Ad (1). The geodesic property follows from the minimality of i3 (cf. with (j0.2|) in 
Definition 14. 4p . Then, there holds 

-^£t(^?(r)) = -p(z?(r), -D£t(i?(r))) = £t(u+) - Et{n-) for a.a. r e (0, 1), (4.10) 
dr 

where the first identity ensues from the chain rule (|2.26[) for £ and the contact condition (|4.3[) . 
and the second one from (|4.2p . Clearly, (|4.10p implies (|4.4p . 

Ad (2). If is of shding type, then the contact condition with ((X^ . yields ((^ : 

follows since d a.e. in (0, 1). 

Ad (3). Equation (|4.7p still follows from (|3.25p . Choosing ro G (0,1) and a Borel selection 
e{r) G A{-d{r), — D£t(?9(r))) (which is therefore locally bounded away from 0), we set 

s(r) := / e-\p)dp, r := s-\ (4.11) 

J 7-0 

so that r is defined in a suitable interval of M and satisfies 

r(.s)-e(r(s)), ^(s) = e(r(s)),9(r(s)). 
Ad (4). We simply introduce the disjoint open sets 

V:=[re (0, 1) : -D£t(79(r)) ^ K*}, S := (0, 1) \V 
and we consider their canonical decomposition in connected components. 

Ad (5). If 7? is energetic, then by dOcncrP and (|4.4p there holds Ap.£ (t; it_, m_|-) = ^'o(u+ — m_). 
Thus, taking into account ([42]) and ([3?3l) as weU, we find p(i?, -D£t(?9(r))) = *o(^5(r)) for a.a. 
r g (0, 1). Since its \l/o-velocity is constant and the total length is \E'o(w-i- — u-), we deduce that 'd 
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is a constant speed minimal geodesic for ^I'o- Conversely, the constraint — D£j(??(r)) G K* satisfied 
by sliding transitions yields, in view of (|3.2f p . that p{'d, — D£j(z?(r))) = 5'o(iS'(r)) for a. a. r e (0, 1). 
Therefore, 

Jo 

by the geodesic property (|4.9p . 

It is well known that, if 4*0 has strictly convex sublevels, the related geodesies are linear seg- 
ments. In order to prove the last statement, let us observe that for every ^ e 9^'o(u+ — u-) C K* 
there holds 

(C, ^{r)) dr = u+ - u_) = ^SJo{u+ - u^) = f *o(i?(r)) dr , 
/o Ja 
where the second equality follows from the characterization (|2.2p of ^'^'o{u-^- — Hence, 



*oWr))-(e,i9(r))jdr-0. 

Since the above integrand is nonnegativc (being ^ G K*), again by (|2.2p we deduce that ^ G 
9*I'o(i?(r)) for a. a. r G (0,1). On the other hand, if '^q is Gateaux-differentiable outside 0, its 
subdifferential contains just one point. Ultimately, (|4.5p (recall that i9 is of sliding type) shows 
that -D£t(i?(r)) = ^ for every r G [0, 1]. □ 

The next result clarifies the relationships between energetic and BV solutions. 
Corollary 4.9 (Energy balance and comparison with energetic solutions). 

(1) A BV solution u of the rate-independent system (X, £,p) satisfies the energy balance (jE| if 
and only if every optimal transition associated with its jump set is of energetic type dOcncrl ) ■ 

(2) A BV solution u is an energetic solution if and only if it satisfies the global stability condition 
([S|. In that case, all of its optimal transition curves are of energetic type. 

(3) Conversely, an energetic solution u is a BV solution if and only if, for every t G J^, any jump 
couple {u{t^),u{t-^)) can be connected by a sliding optimal transition. 

Proof. Ad (1). Let m be a BV solution such that every optimal transition is of energetic 
type dOcncrl )- Now, taking into account ( |Jbv[ ), one sees that dOcncrl ) is equivalent to the jump 
conditions dJcnciD - Then, equation ( DNq.bvP (which holds by Theorem 14. 3p and dJcncrD yield the 



energy balance ([E| (cf. the proofs of Propositions 12.21 and 12. 7p . The converse implication ensues 
by analogous arguments. 

Ad (2). The necessity is obvious; for the sufficiency we observe that, for every jump point t G Ju, 
the global stability condition jS]) (written first for it(t_) with test functions v = u(t) and v — uit-^-), 
and then for uit) with v — u{t^)), yields 

*o(u(i) -M(i-)) > £t(w(t-))-£t(M(t))= Ap,£(i;u(i_),u(i)) > ^>oiuit) - u{t,)), 
Mu{t+)-u{t.))> £t{u{t-)) - £t{u{t+)) = Ap,£(t;w(i_),u(t+)) > ■^o{u{t+) - u{t.)), 
'fo{u{t+)-u{t))> £t{u{t)) - £t{u{t+)) ^ Af,x{t;u{t),u{t+))> yfo{u{t+)~u{t)), 

where the intermediate equalities are due to dO-2p and the subsequent inequalities to (|3.9[) . The 
resulting identities ultimately show that the transition is energetic, by the very definition dOcncrD - 
Ad (3). The condition is clearly sufficient. It is also necessary by the previous point, since 
energetic transitions are in particular of sliding type. □ 

4.3. Viscous limit. We conclude this section by our main asymptotic results: 

Theorem 4.10 (Convergence of viscous approximations to BV solutions). Consider a sequence 
(we) C AC([0,T];X) of solutions of the viscous equation l \DN^l , with u^{0) — > uq as £ J, 0. 

Then, every vanishing sequence £fc | admits a further subsequence (still denoted by (sk) ), and 
a limit function u G BV([0,T];X) such that 

Uef,(t) — > u{t) for every t G [0, T] as k 1 +oo, (4.12) 
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and u is a BV solution of ( |DNo[ ), induced by the vanishing viscosity contact potential p according 
to Definition \4-.l\ 



Proof. It follows from the discussion developed in Section 12.51 that for every sequence Sk [ there 
exists a not relabeled subsequence [u^^,) such that (|4.12[) holds, and u complies with the local 
stability condition dSiod ). In view of Proposition 14.21 it is then sufhcient to check t hat (|Ep £:iiieq ) 
holds. The latter energy inequality is a direct consequence of the e-energy identity (|2.27p and the 
lower semicontinuity property stated in Lemma 16.151 later on. □ 

Our next result concerns the convergence of the discrete solutions to the viscous time-incremental 
problem ( |IPg[ ), as both the viscosity parameter e and the time-step r tend to 0. 



Theorem 4.11 (Convergence of discrete solutions of the viscous incremental problems). Let 
Ut.e : [0,T] X be the left- continuous piecewise constant interpolants of the discrete solutions of 
the viscous incremental problem ( |IPe[ ), with U° ^ — > uq as e, r | 0. 
Then, all vanishing sequences Tk , £fc J. satisfying 

lim— = +oo (4.13) 

admit further subsequences (still denoted by (t^) and (sk)) and a limit function u G BV([0, T]; X) 
such that 

Urfc,efc(i) u{t) for every t e [0, T] as k ] +oo, 

and u is a BV solution of ( |DNo[ ) induced by the vanishing viscosity contact potential p according 
to Definition\^.l\ 



The reader may compare this result to [TBI Ell [22J [49] , where the same double passage to the limit 
was performed for specific applied problems and conditions analogous to ()4.13p were imposed. 



Proof. The standard energy estimate associated with the variational problem (jlPJ) yields 

-^^[^{Ul\, - f/;-^)) + £-t,AK\,) < £-tAK-') = £*„_,([/;-!) + £" dt£t{Ke')dt. (4.14) 

Thanks to (|2.7p , we easily get from (|4.14[) the following uniform bounds for every 1 < n < (here 
C is a constant independent of n, r, e) 

N N 



e \T 

n—l 



the latter estimate thanks to (|^.2|) . 

Denoting by ^ (resp. Ur.e) the right- continuous piecewise constant interpolants (resp. piece- 
wise linear interpolant) of the discrete values {U"^) which take the value [/"^ at t = tn, we have 

£t(Ur,e(t)) < C, Var>t„(U,,,; [0,T]) < C (4.15a) 
||U,,, -lJ,,e||L-(o,T;X), ||U,,e "U.,J|l="(o,T;X),< supllt/;;, - K:;'\\x < Cujir/Ce), (4.15b) 

n 

where 

uj{r) := sup {\\x\\x : r^{r^'^x) < l} 

satisfies limrxo'^('') = thanks to (|^.ip . By Holly's theorem, these bounds show that (up to 
the extraction of suitable subsequences (rk) and (e^) satisfying (|4.13p ). the sequences (Urk,ek)j 
(Ui-fc,efc) and (Htj. gj.) pointwise converge to the same limit u. 

By differentiating the variational characterization of C/"^ given by ( |IPg[ ) we obtain 

d^J^ 1:^) + 3 0, -D£,„ ([/"), 
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which yields in each interval (tn-ii^n] (here, Wr^e denotes the left-continuous piecewise constant 
intcrpolant of the values {W"^)^^i) 

where 

R{t; X, y) ~ £t(y) - lt{x) - {D£t{y),y - x). 
Since £ is of class C^, for every convex and bounded set B <Z X there exists a concave modulus of 
continuity ub ■ [0, +00) [0, +00) such that linirjo <7B{f) ~ crs(O) = and 

R{t\x,y) < (Tb{\\v - x\\x)\\y - x\\x for every t e [0,r], x,y e B. 

We thus obtain 

i-T 



(^',(U,,aO) + *:(W,,eW))di + £t„(U,,e(iiv)) < £oM + J 9t£t(U,,,(t))dt 

N 

sup asiWK., - ^ ~ WrAt)) = -D£t^(,)(lJ.,,(i)). 



(4.16) 



We pass to the limit along suitable subsequences (rfe) and {sk) such that Ur^.c^, Ur^.e^ u 
pointwise; since Ut,£ and IJt,s are uniformly bounded, (|4.15bp and ()4.13p yield the convergence to 
of the third term on the right-hand side of (|4.16p . which thus tends to £0(1*0) + Jq 9t£t(u(t)) dt. 
Since Wrfc^£j.(i) — > w{t) = — D£t(u(t)), applying the lower semicontinuity result of Lemma [6.151 



we obtain that u satisfies (Ep.£:incql and the local stability condition. In view of Proposition 14.21 



this concludes the proof. □ 

5. Parametrized solutions 



In this section, we restart from the discussions in Sections 12.41 and 12.51 and adopt a different 
point of view, which relies on the rate-independent structure of the limit problem. The main idea, 
which was introduced by [18j , is to rescale time in order to gain a uniform Lipschitz bound on the 
(rescaled) viscous approximations. Keeping track of the asymptotic behavior of time rescalings, 
one can retrieve the BV limit we analyzed in Section 3) In particular, we shall recover that the 
limiting jump pathes reflect the viscous approximation. 

5.1. Vanishing viscosity analysis: a rescaling argument. Let us recall that for every e > 
Me are the solutions of the viscous differential inclusion 

a«'e(ite(i)) + D£t(ue(t)) 3 in X* for a.a. t G (0, T), (DN^) 

which we split into the system 

D£t(Ue(t)) = -We, dt£.t{Ue{t)) = -Pe- 

We follow the ideas of [18l [35j to capture the aforementioned limiting viscous jump pathes,. 
However, owing to the dissipation bound (|2.28p . we use a different time rescaling Se : [0, T] [0, Se] 

Seit) + (*c(we(r)) + «':(iOe(r))) dr and Se := Se(T). (5.1) 

Thus, Se may be interpreted as some sort of "energy arclength" of the curve Me. Notice that, 
thanks to (|2.28p . the sequence (Se) is uniformly bounded with respect to the parameter s. Let us 
consider the rescaled functions (tg, Ue) : [0, Se] — > [0, T] x X and (pe, We) : [0, Se] ^ M x X* defined 
by 

te(s) :=S-l(s), Ue(s) :=Me(te(s)), 

(5.2) 

Pe(s) :=Pe(te(s)) = -5t £t^ (Ue (s)) , We(s) :=M;e(te(s)) = -D£t^(s) (Ue(s)). 
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We now study the limiting behavior as e 4 of the reparametrized trajectories 

{(te(s),u,(s)) : s e [0,Se]} C = [0,r] X X, 
{(te(s),ue(s);pe(s),w,(s)) : se[0,Se]}c^, 
where we use the notation 

^ := [0,+oo) X X X M X X*. (5.3) 

In order to rewrite the "rescaled energy identity" fulfilled by the triple (te, UejW^), we define the 
viscous space-time vanishing viscosity contact potential : (0, +oo) x X x M x X* [0, +oo) 
by setting 

as ot 

'^^(a, v; p, w) := ^^'^(v/a) + a^'*(w) + ap = -'I'(-v) + -'I'*(w) + ap (5.4) 

e a £ 

Hence, (|2.27p becomes for all < si < S2 < 

S2 

*Pe(t£(s),Ue(s); Pe(s),W£(s))ds + £t^(,^) (Ue(s2)) = £t^ (^^) (u,. (si)), (5.5) 

L 

and (|5.ip yields 

q3e(te(s),Ue(s);l,We(s)) =1 for a.a. s G {G,Se) ■ 

A priori estimates and passage to the limit. Due to estimate (|2.28p . there exists S > such 
that, along a (not relabeled) subsequence, wc have Se{T) S as e J, 0. Exploiting again (|2.28p . 
the Arzela-Ascoli compactness theorem, and the fact that X is finite-dimensional (see also the 
proof of [351 Thm. 3.3]), we find two curves t G W^'°°{0,S) and u G W^'°°{[0,S];X) such that, 
along the same subsequence, 

te ^t in C"([0,S]), te^*t in L°°(0,S), (5.6a) 

Ue ^ u in C"([0,S];X), u^^* u in i°°(0, S; X), (5.6b) 

Pe ^ p in C°([0,S]), We ^ w in C°([0,S];X*), (5.6c) 

with 

£t.(ue)^£t(u), p(.s) = -9t£t(,,)(u(.s)), w(.s) = -D£t(,)(u(s)) (5.6d) 
for all s G [0, S]. Then, to pass to the limit in (j5.5p we exploit a lower semicontinuity result 
(see Proposition 16. 2p . based on the fact that the sequence of functionals (^e) F-converges to the 
augmented vanishing viscosity contact potential Cp : [0, -f oo) x X x M x X* — > [0, +oo] (see Lemma 
I6.ip defined by 

X /*o(«)+lK-(w)+ap ifa>0, 
^{a,v;p,w) i (5.7) 
\p(v, w) if a = 0. 

By (|5.6p and Proposition 16. 21 we take the liminf as e J, of (|5.5p and conclude that the pair (t, u) 
fulfils, for all < si < S2 < S, the estimate 

^{i{s), u{s); p(s), w(,s)) d,s + £t(.,)(u(s2)) < £t(.,)(u(si)) • (5-8) 

5.2. Vanishing viscosity contact potentials and rate-independent evolution. The aug- 
mented space-time contact potential *p is closely related to p introduced by (|3.13p . The following 
result fixes some properties of Its proof, which we choose to omit, can be easily developed 
starting from Theorems 13.61 and 13 . 121 for the vanishing viscosity contact potential p. 

Lemma 5.1 (General properties of ?P). 

(1) ^ is lower semicontinuous, 1-homogeneous and convex in the pair (a, f); for every {ol,v) G 
[0, -|-oo) X X the function Cp(a, w; ■, ■) has convex sublevels. 

(2) For all {a,v,p,w) G ^ (cf. (|5.3p ) it satisfies 

'^{a,v;p,w) > {w,v) + ap, *P(0, w) > p(?;, w) > Vl'o(w), (5.9) 
^{Q,v]p,w) ^-^oiv) ^ weK*. (5.10) 



26 



ALEXANDER MIELKE, RICCARDA ROSSI, AND GIUSEPPE SAVARE 



(3) The contact set of *p 

S(p := |(q!, v;p,w) € ^ : *P(q!, v;p, w) = {w, v) + ap| (5-11) 
does not impose any constraint on p. It can be characterized by 

(q;,-u;p, w) e Sqj ^ w G d^{a, ■ ]p,w){v). (5-12) 

We also have 

for a > 0, {a,v;p,w) G Stp if and only if w £ 9^'o(w), (5.13) 

for a ~ 0, (a, w) G Eq3 if and only if (v,w)£'Sp. (5-14) 
Equivalently, {a,v;p,w) G Stp if and only if 

w e d'i'oiv) C K* or [w^K*, a = 0, 3e E A{v,w) : we^^^{ev)], (5.15) 



where A{v,w) is defined in (|3.18p . In particular, in the additive viscosity case (j2.2ip . we 
simply have 

{a,v;p,w) eY^<^ 3A>0: w G a^'o(w) + S^fy (Aw) and aA = 0. (5.16) 

Conclusion of the vanishing viscosity analysis. We are now going to show that (|5.8p is 
in fact an equality. This can be easily checked relying on the chain rule (|2.26p . which yields 
for a. a. s G (0, S) 

;^£t(,)(u(s)) = -a,£t(.)(u(s))t(s)-(-D£t(.)(u(s)),u(.)) 

ds ' ' (5.17) 

^ ^p{s)i{s) - (w(s), u(s)) > -^{{{s), u(s), p(s), w(s)) . 

Collecting (|5.17p and (|5.8p . we conclude that the latter holds with an equality sign and. with an 
elementary argument, that such equality also holds in the differential form, namely for a. a. s G 
(0,S) 

p(s) = -at£t(s)(u(s)), w(s) = -D£t(s)(u(s)) 
^£t(,)(u(s)) = -p(s)t(s) - (w(s),u(s)) = -<P(t(s),u(s),p(s),w(s)) ^^■^^^ 

which yields 

(t(s), u(s); ~at£t(,) (u(s)), -D£t(,) (u(s))) G S<p for a.a. s G (0, S). (5.19) 
Finally, we take the lim sup as e | of (|5.5p , using (|5.6p and (|5.18p , whence 

limsup / <Pe(te(s),Ue(s),Pe(s),We(s))ds < / «P(t(s),u(s),p(s),w(s))ds. 

eiO Jo Jo 

In particular, we find that for a.a. s G (0, S) 

<P(t(s),u(s);l,w(s)) = 1. (5.20) 

5.3. Parametrized solutions of rate-independent systems. Motivated by the discussion of 
the previous section, we now give the notion of parametrized rate-independent evolution, driven 
by a general vanishing viscosity contact potential satisfying conditions (1), (2) of Lemma [5. II 

Definition 5.2 (Parametrized solutions of rate-independent systems). Let *p : ^ — > (—00,-1-00] 
be the vanishing viscosity contact potential (j5.7p . We say that a Lipschitz continuous curve 
(t, u) : [a,b] — > [0,T] x X is a parametrized rate-independent solution for the system (X, £,*P) if 
t is nondecrcasing and, setting p(s) = —dtEt{s){^{s)), w(s) = — D£t(s') (u(s)) for all s G [a,b], we 
have 

q3(t(s), u(s); p(s), w(s))ds + £t(3,) (u(s2)) < £t(si) (u(si)) Va < si < S2 < b. (5.21) 



Furthermore, 

(1) if t(s) -I- ^'o(u(s)) > for a.a. s G (a, b) we say that (t, u) is nondegenerate; 
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(2) if t(a) = 0,t(6) = T we say that (t, u) is surjective; 

(3) if (t, u) satisfies (j5.20p . we say that it is normalized. 



Definition 15.21 generalizes to the present setting the notion whieh we first introdueed in [35| . 

Remark 5.3. The nice feature of the previous definition is its invarianee with respect to (nonde- 
creasing, Lipschitz) time rescalings. Namely, if (t, u) : [a, b] — > [0, T] x X is a parametrized solution 
and s : [a, f3] — >• [a, b] is a Lipschitz nondecreasing map, then (t o s, u o s) is a parametrized solution 
in [a,P]. 

The next result provides equivalent characterizations of parametrized solutions. 

Proposition 5.4. A Lipschitz continuous curve (t, u) : [a,b] — > [0, T] x X, with t nondecreasing, 
is a parametrized solution of {X, £,^P) if and only if one of the following (equivalent) conditions 
(involving as usual p = — 9t£t(u),w = — D£t(u)j is satisfied: 
(1) The energy inequality (|5.2ip holds just for si = a and S2 = b, i.e. 

rb 

^P(t(s), u(s); p(s), w(s))ds + £t(b)(u(6)) < £t(a)(u(a)). (5.22) 



(2) The energy inequality ()5.2ip holds in the differential form 

^£t(,)(u(s)) +<p(t(s),u(s);p(s),w(s)) < for a.a. s € ia,b). (5.23) 

(3) The energy identity holds, in the differential form 

-^£t(,)(u(.s))+*P(t(s),u(s);p(5),w(s)) = for a.a. s e {a,b), (5.24) 
as ^ ' 

or in the integrated form 

<P(t(s),u(s);p(s),w(s))ds + £t(32)(u(s2)) = £t(si)(u(si)) for a < si < S2 < b. (5.25) 

(4) There holds 

(t(s),u(.s);-at£t(s)(u(s)),-D£t(s)(u(s))) e Sip for a.a. s <E {a,b) . 

(5) The pair (t, u) satisfy the differential inclusion 

a'P(t(s),.;"5t£t(,)(u(s)),-D£t(,)(u(s)))(u(s))+D£t(,)(u(s)) 9 a.e. m {a,b) . (5.26) 

In particular, for a.a. s £ {a,b) we have the implications 

t(s)>0 => -D£t(,)(u(s)) e A-, 
-D£t(,)(u(s)) e K* => -D£t(,)(u(s)) e a*o(u(s)), 

and for every Borel map A defined in the open set J by 

J:={,se(a,5):-D£t(,)(u(s))^A-*}, 

with A(s) G A(u(s), — D£t(s) (u(.s))) for a.a. s G J, 

we have 

-D£t(,)(u(s)) e 5^'(A(s)u(s)), t(s) = for a.a. s e ]. (5.29) 

The proof follows from the chain rule ^(TI^ (arguing as for ([STT]) . ((57T5)) . ([STTO]) ). and from 
the characterization of the contact set Sqj of Lemma [5TT] (see also [3S1 Prop. 2]). 

Corollary 5.5 (Differential characterization in the additive viscosity case). Let ^ : S§ —> (— oo, +oo] 
be a vanishing viscosity contact potential satisfying conditions (1), (2) of Lemma \5.1[ and suppose 
also that the contact set of *P satisfies the characterization (|5.16[) of Lemma \5.1\ in the additive 
viscosity case (|2.2ip * = *o + ^'v- 

Then, a Lipschitz continuous curve (t, u) : [a, 6] — > [0,T] X is a parametrized solution of 
{X, £, *P) if and only if there exists a Borel function A : (a, b) [0, +oo) such that for a.a. s G (a, b) 

a*o(u(s)) +a*y(A(s)u(s)) + D£t(,)(u(s)) 9 0, A(s)t(s) = /or a.a. s e (a, 6). (5.30) 



(5.27) 



(5.28) 
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The vanishing viscosity analysis developed in Sections 15.11 and 15.21 provides the following con- 
vergence result. 

Theorem 5.6 (Convergence to parametrized solutions). Let (u„) be viscous solutions of ( |DNg[ ) 
corresponding to a vanishing sequence (£«), let t„ : [0,S] —>■ [0,T] be uniformly Lipschitz and 
surjective time rescalings and let u„ : [0, S] X be defined as u„(s) := Un(tn(s)) for all s € [0,S]. 
Suppose that 

3a >0 VneN : m„(.s) := (t„(s), u„(s); 1, -D£t„(,) (u„(s))) G [a, a"!] 

for a.a. s S (0, S). If the functions (tn, u„, m„) pointwise converge to (t, u, m) as n oo, then (t, u) 
is a (nondegenerate, surjective) parametrized rate-independent solution according to Definition \5.'A 
and 

^5(t(s), u(s); 1, -D£t(,) (u(s))) = m(s) for a.a. s e (0, S). 

The following remark, to be compared with Remark 14.81 highlights the different mechanical 
regimes encompassed in the notion of parametrized rate-independent solution. 

Remark 5.7 (Mechanical interpretation). The evolution described by (|5.26p in Proposition 15.41 
bears the following mechanical interpretation (cf. with [18] and [35]): 

• the regime (t > 0, ii = 0) corresponds to sticking, 

• the regime (t > 0, u 7^ 0) corresponds to rate-independent .sliding. In both these two 
regimes -D£t(u) £ K* . 

• when — D£t(u) cannot obey the constraint /\*, then the system switches to a viscous 
regime. The time is frozen (i.e., t = 0), and the solution follows a viscous path. In the 
additive viscosity case (|2.2ip it is governed by the rescaled viscous equation (|5.30p with 
A > 0. These viscous motions can be seen as a jump in the (slow) external time scale. 

We conclude this section with the main equivalence result between parametrized and BV so- 
lutions of rate- independent systems (compare with the analogous [351 Prop. 6]). We postpone its 
proof at the end of the next section. 

Theorem 5.8 (Equivalence between BV and parametrized solutions). Let (t, u) : [0, S] [0, T] x 

X be a (nondegenerate, surjective) parametrized solution of the rate independent system (X, £,*P). 
For every t € [0, T] set 

s(i) := {s G [0,S] : t(s) = t} (5.31) 
Then, any curve u : [0,T] X such that 

u(t) G {u(s) : s G s(i)} (5.32) 
is a BV solution of the rate-independent system {X, £,p). 

Conversely, if u : [0,T] X is a BV solution, then there exists a parametrized solution (t, u) such 
that (|5.32p holds for a time-res caling function s defined as in (|5.3ip . 

6. Auxiliary results 

After proving some lower semicontinuity results for vanishing viscosity contact potentials, in 
Section 16.21 we develop some auxiliary results concerning the total variation induced by time- 
dependent (and possibly asymmetric) Finsler norms. 

6.1. Lower semicontinuity for vanishing viscosity contact potentials. Let us start with a 
lemma which shows that which is defined in (|5.4p . F-converges to *P as e 4 (compare with 
[35l Lemma 3.1]), where *p is defined in (|5.7p . 

Lemma 6.1 (F-convergencc of ^e). 

F-Hminf estimate: For every choice of sequences £„ ], and (««, fnjPnj w„) — > {a,v,p,w) 
in SS , we have 

liminf <P£„(a„,w„;p„, w„) > ^i{a,v;p,w). (6.1) 
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F-limsup estimate: For every {a,v;p^w) G there exists (a^, u'e)e>o such that 

limsup*Pe(ae,We;pe, We) <^{a,v;p,w). (6.2) 

eiO 

Proof. The F-liminf estimate is easy: if a > then, also recahing (|2.29p . one verifies that 
hminf <P£^(a„,z;„;p„, Wn) > liminf (*o('^«) + ane:^'^%{wn) + a„p„) 

n— >oo n— »oo ^ (6 3) 

> ^o{a,v;p,w), 

where we have used the notation 

<Po(a, v; p, w) := *o(«) + Ik{w) + ap. (6.4) 

The first inequaUty in (|6.3p is also due to (|2.f 6p . If a = 0, we use the obvious lower bound 

*Pe„(a„,w„;p„,-u;„) > p(w„,-u;„) + a„p„ 

and the continuity of p (cf. Theorem 13. 12|) . 

To show the limsup estimate (j6.2p for w G K*, we simply choose := a + e, ;= v,p^ := 
p, We := w, observing that in this case 

*Pe(ae, We) < e{a + e)-^{v/{e{a + e)) + {a + e)p "-^ ^'o(w) + ap = *P(a, w) , 

the first passage due to (|2.24p . If w ^ K* , we choose a coefficient A G A(w, w) as in ()3.I8p . and we 
set := Ae, := v,Pe := p, := w, obtaining 

VE{o!E,Ve;PE,We) = p{v,w) + Xep "-^ p{v,w) = <P(a,t;;p, w). □ 

An important consequence of the previous Lemma is provided by the following lower-semicontinuity 
result for the integral functional associated with 

Proposition 6.2 (Lower-semicontinuity of the e-energy). Let us fix an interval (sq, si). For every 
choice of a vanishing sequence > and of functions a„ G i°°(so,si), p„ G L^{so, si), w„ G 
L'^{0,T;X), Wn e L^{0,T;X*) such that 

a„— ^* a in L°° {sq, si), Pn ^ P in L^{0,T), 
Vn ^ V in L^{0,T;X), Wn ^ w inL^{so,si), 

we have the liminf estimates 

liminf / ^Pe„ (a„(s), w„(s);p„(s), w„(s)) ds > / ^{a{s),v{s);p{s),w{s)) ds, (6.5) 

J So J So 

rsi rsi 

liminf / *Po(ari(s), u„(s);p„(s), w„(s)) ds > / ^{a{s),v{s);p{s),w{s)) ds, (6.6) 

J so J So 



where defined in (|6.4p . 

Proof. It is sufficient to prove this result in the case pn = p — Q. Then we notice that, by Lemma 
16.11 the integrand 

*P(e, a, V, w) := v; 0, w) for (e, a, w, w) G [0, +oo) x [0, +oo) x X x X* 

is lower semicontinuous and convex in the pair (a,w). Then, inequality (j6.5p follows from loffe's 
Theorem (see e.g. [31 Thm. 5.8]). A similar argument yields (|6.6p . □ 

6.2. Asymmetric dissipations, pseudo-total variation, and extended space-time curves. 
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Notation. Hereafter, St^ shall stand for the extended space-time domain [0, T\ x X, with elements 
X = (t, u) denoted by bold letters. We shall denote by "V the tangent cone [0, +C!o) x X to ^ and 
V — (a, v) the elements in . 

We shall consider lower semicontinuous dissipation functionals 3?: ^ xY ^ [0,+cxd) satisfying 
the following properties: 

Va; G ^ : 3? (a;; •) is convex and positively 1-homogeneous; (6.7a) 
3C> Va; G JT, = e ^ : ■Jl{x-v)>C\\v\\x (6.7b) 

is lower semicontinuous on ^ x Y . (6.7c) 

In order to keep track of the time-component of v we also set, for all /3 > 0, 

3?^j {x\v) ^ aj3 + ']l{x\v) for aUa; G JT, v = {a,v) G Y. 

Notice that, for any dissipation 'Ji complying with properties (j6.7[) . the corresponding functional 
'Jlp satisfies the subadditivity property for &\\ x ^ and Vi, V2 £ Y 

Jip {x; Vi + V2) < Rp {x; Vi) + 3?^ {x; V2) . 

Example 6.3 (Dissipations induced by ^Pq and 

(1) Our first trivial example of a dissipation fulfilling properties ()6.7p is given by 

fix, v) := l'o(u) for a; G jr, v = (a, v) G t. (6.8) 

(2) Our main example will be provided by the dissipation induced by the vanishing viscosity 
contact potential namely 

'B{x;v) ■.= f{a,v,Q,-T)E,t{u)) for x ^ {t,u) E ^ , v ^ {a,v) & Y . (6.9) 



It is not difficult to check that 23 satisfies all of assumptions (|6.7p . Hence, for all /? > we 

set 

"Bisix-v) ■.= ^{a,v;(3,-J:)£.t{u)) for x ^ {t,u) e ^' , v = {a,v) G Y . (6.10) 

Definition 6.4 (Pseudo-Finsler distance induced by D?). Given a dissipation function 51 : ^ x 
Y [0, +00) complying with ([ej]) . for every xi = {U,Ui) G i = 0,1, with < to < < T, 
we set 



Aa^^ (a;o, a;i) := inf / Jlp (x(r); x(r)) dr : 

^-''■0 (6.11) 

X = (t, u) G Lip(ro, ri; JT), a?(ri) = a;,;, i = 0, 1, t > o|. 
If to > ti we set A,'Kg(a;o,a;i) := -l-oo. We also define 

Aoio{t;uo,ui,U2) := Aoio{{t,uo), {t,ui)) + AKo((t,wi), {t,U2)). 
(notice that this quantity is independent of /3). 

Remark 6.5. The link with the Finsler cost Ap_£ (|3.8p induced by (p, £) is clear. For H = H given 
by (|6.9p . using *P(0, w; 0, w) = p(?;, w) we have, for to = ti = t, 

A'Bo{{t,uo),{t,ui)) = Apxit;uo,ui) for every uo,ui G X. (6-12) 

When 3? = CP is given by (|6.8p . we simply have 

Ay^,{{to,uo),{ti,ui)) = /3(ti - to) + «'o(wi - Wo) for uo,ui e X, < to < < 
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General properties of A^^{-, •). It is not difRcult to check that the infimum in (|6.1ip is attained 
and, by the usual rcscaling argument (cf. Remark l4.5p . one can always choose an optimal Lipschitz 
curve X = (t, u) defined in [0, 1] such that 

3^1 (x; x) is essentially constant and equal to A^^ {xo,xi) ~ (ti — to) + A^r^ {xq, xi). (6.13) 

Properties (j6.7bp - (|6.7cp yield, for every uo,ui G X and < < ti < T, the estimate 

Pih - to) + C\\ui - uoWx < A3j^((io,uo),(ii,ui)). (6.14) 

Notice that Aji.^{-, ■) is not symmetric but still satisfies the triangle inequality: for Xi = {ti,Ui) e 
^ with to < ti < t2, there holds 

Aoi^{xo,X2) < Aoii^{xo,xi) + Aoif^{xi,X2). 

Another useful property, direct consequence of (|6.7c[) . is the lower semicontinuity with respect to 
convergence in if Xi_n = (ti^„ , Ui^n) Xi = {ti, Ui) in j?r as n f +oo, i = 0, 1, then 

liminf Ajj (a;o,„,a;i,„) > Ajj^ (a;o, a;i). (6.15) 

Indeed, assuming that the liminf in (|6.15p is finite and that, up to the extraction of a suitable 
subsequence, that it is a limit, it is sufficient to choose an optimal sequence x„ = (t„,u„) of 
Lipschitz curves as in (|6.13p . which therefore satisfies a uniform Lipschitz bound and, up to the 
extraction of a further subsequence, converges to some Lipschitz curve x = (t, u). Then, (|6.15p 
can be proved in the same way as (|6.6p . 

In the case of the cost induced by 23 induced by the vanishing viscosity contact potential *P, we 
have a refined lower-semicontinuity result: 

Lemma 6.6. Let Un,Wn ■ [^07^^11 X be Borel maps and (e„) be a vanishing sequence. Suppose 
that u.fi is absolutely continuous for every n € N and that the following convergences hold as n oo 

Un{t) ^ u(t) and Wn(t) w(t) for all t ^ [to,ti]; sup ||w„(i) + D£t(u„(t))||x' — * 0. 

te[to,ti] 

Then, 

(^^,^iM-t)) + K.M))) dt > Aji,iito,uito)),it,,u{ti))). (6.16) 



lim inf 

)iT+oo Jf.^ 



Proof. Up to extracting a further subsequence, it is not restrictive to assume that the liminf in 
([Oe]) is in fact a hmit. We set as in jSH]), 

Sn{t) -.^ t - to + (*e„ «(?■)) + *e„(«^"(0)) dr, S„:=s„(ti), 
t„(s) := s,7^(s), u„(s) := u„(t„(s)), w„(s) := w„(t„(s)) for all s G [0, S„] 

so that 

(*,„«(i)) + *:„(«i„(t)))dt = ^ "'P,„(t„(s),u„(s);0,w„(s))ds. (6.17) 

Since the sequences (t„) and (u„) are uniformly Lipschitz, applying the Ascoli-Arzela Theorem we 
can extract a (not relabeled) subsequence such that S„ — > S, and find functions t : [0, S] —>■ [to, ti], 
u : [0, S] X, and w : [0, S]-^ X*, such that 

t„ t, u„ u, w„ ^ w = — D£t(u) uniformly in [0, S]. 

By construction, we have t(0) = to, u(0) = u{to), t(S) = ti, and u(S) = u{ti). Applying Proposi- 
tion [621 we have 



lirninf / <Pe„ (t„(s), u„(s); 0, w„(s)) ds > / <Po(t(.s), u(.s); 0, w(s)) ds 

"T+oo Jg Jo (6.18) 

> A3o((t(0),u(0)),(t(S),u(S))). 
Combining (|6.17p and (|6.18p . we conclude (|6.16p . □ 
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The total variation associated with Aj^.^. In the same way as in Definition 13.41 we introduced 
the total variation Varp^g associated with the Finsler cost Ap g, it is now natm'al to define the 
total variation associated with Agj^, . 

Definition 6.7 (Total variation for the pseudo-Finsler distance Ajj^). For every curve x = (t, u) : 
[0, S] ^ ^ such that t is nondecreasing and every interval [a, b] C [0, S] we set 

M 

Var5;^(x; [a, b]) :=sup| Ajj^^ (x(s„J, x(s™_i)) : 

^ ™=i (6.19) 

a = So < si < • • • < SA/-1 < SM = b|. 

For a non-parametrized curve u : [0,T] X and [a, b] C [0, T], we simply set 

Var3j^(w; [a,5]) := VarK^(M; [a,6]), with u{t) {t,u{t)) e ^ , t e [0,T]. 

In view of (j6.7bp . it is immediate to check that a curve u with Ya.r^g{u; [0, T]) < +oo belongs 
to BV([0,T];X). 

In contrast to the (pseudo)-total variation defined in p. lip , the above notion of total variation 
is lower semicontinuous with respect to pointwise convergence (compare with Remark 13. 5|) . 

Proposition 6.8 (Lower semicontinuity of Vargj^, (•; [a, b])). // x„ = (t„, u^) : [0, S] ^ ^ is a 

sequence of curves pointwise converging to x ~ (t, u) as n ^ oo, we have 

liminf Var3j,(x„; [a, b]) > Var3j„(x; [a, b]). (6.20) 

nfoo 

Proof. The argument is standard: for an arbitrary subdivision a = sq < si < • ■ • < sm-i < sm = 
b, ([635]) yields 

M M 

> A3j^(x(sm),x(sm_i)) < liminf > A3j^(x )) < liminf VarK^(x„; [a,b]). 

^ — ^ nj+oc ^ — ^ ntoc 

m— 1 m— 1 

Taking the supremum with respect to all subdivisions of [a, b] we obtain (|6.20|) . □ 

Lipschitz curves. The next result shows that, for Lipschitz curves, the total variation can be 
calculated by integrating the corresponding dissipation potential. 

Proposition 6.9 (The total variation for Lipschitz curves). Given (3, L > 0, a bounded curve 
X := (t, u) : [0, S] —^ ^ satisfies the Ajj^ -Lipschitz condition with Lipschitz constant L 

Ak^(x(si),x(s2)) < L(s2 - si) for every < si < S2 < S, (6.21) 

if and only if it is Lipschitz continuous (with respect to the usual distance in ^ ), t is nondecreasing, 
and 

Jlp (x(s); x(s)) < L for a.a. s G (0, S). (6.22) 

In this case, for every 7 > 

Var3j^(x;[a,b])=7(t(b)-t(a))+^ 3^0 (x(s); x(s)) ds. (6.23) 

Proof. The sufficiency of condition (|6.22p is clear. Let us now consider a curve x satisfying (|6.21[) : 
by the coercivity (|6.14p . a; is a Lipschitz curve in the usual sense and [47l Prop. 2.2] yields 

Aoi^{x{s),xis + h)) 



Ajj , (x(.so), < / m(s)ds, where m(s):=lim- 



h 

is the so-called metric derivative of x (see [US]). The minimality of m ensures that 

m(s) < (x(.s); x(s)) for a.a. s G [0, S]. (6.24) 

On the other hand, since Jip is lower semicontinuous and 1-homogeneous in v, for every < a < 1 
and s e [0, S] we find a constant 6 > such that 

5if3{x{r);v)>a:Rp{x{s);v) ior every v & f if\r-s\<S, 
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SO that a comparison with the Uncar segment joining x(s) with x(s + h) yields 

A3j,(x(s),x(s + /i)) < a-^Jlfi{x{s);x{s + h)-x{s)) 

Dividing by h and passing to the hmit first as ft- 4 and eventuaUy as cr | 1, we obtain the opposite 
inequahty of (|6.24p . Combining (|6.24p (which holds as an equality) with (|6.2ip . we infer (|6.22p . 
and (|6.23p ensues. □ 

Proposition 6.10 (Reparametrization). Let u : [0,T] X be a curve with finite total variation 
V :— VarKo(w; [0, T]) < +oo, and let us set 

s{t) := t + VarKo(u; [0,i]) = VarKi(w; [0,t]) for every t G [0,T]. (6.25) 

Then, there exists a Lipschitz parametrization x = (t, u) : [0, S] 3^ , with S = V + T , such that 

3li (x(s);x(s)) = 1 for a.a.s E {0,S), (6.26) 

t{s{t)) = t, u{s{t)) = u{t) for every t e [0,T]. (6.27) 

In particular, 

Mb) 

b-a + YaToi„{u;[a,b]) = s{b)-s{a)= 3?i (x(s); x(s)) ds. (6.28) 

Js{a) 

Proof. The proof is classical, at least when the dissipation 51 is continuous and even in its second 
argument: we briefly sketch the main ideas and refer to [35|, Lemma 4.1]. 

Notice that the jump set Js of the curve s given by ()6.25p coincides with the jump set J„ of u, 
and s is injective in Cu := (0,r) \ Ju- We denote by t its inverse, defined on := s(C„) and 
extended to Cu by its (Lipschitz) continuity; we also set u(s) := u{t) if s = s{t) e Cu. Suppose 
now that (s_,s+) is a connected component of [0, S] \ Cu, corresponding to some time t g [0,T] 
with s± = s{t±) and s ~ s{t) G [s_, s+J. We have 

u(s_) = lim u(s) = u(t_), u(s+) = lim u(s) = u{t+), 

s - s_ = A3j„ ((?,«(<_)), it,u{t))), s+-s = A3j„((i, «(<)), it,uit+))). 

By Definition 16.41 we can join {t, u(s_)) to {t, u(s+)) by a Ajj^ -Lipschitz curve (still denoted by 
(t, u)) defined in [s_,s_|-] with constant first component t(s) = t, and satisfying (|6.13p as well as 
u(s) = u{t). 

It is then easy to check that the final curve x = (t, u) obtained by "filling" in this way all the (at 
most countable) holes in [0, S] \ C„ satisfies (j6.27p and the Lipschitz condition (|6.2ip with L < 1. 
Applying and ([g?^ we get 

ib) 

Jil (x(.s); x(s)) ds < s{b) - s{a) = YaToi,{u; [a, b]) < Var^i (x; [s(a), s(6)]) 

(a) 

(fc) 

5li (x(s);x(s)) ds 

5(a) 

where the first inequality follows from the 1-Lipschitz condition, the subsequent identity from the 
definition of s, and the last one from (|6.28p . □ 

The reparametrization of Proposition lS.lOl is also useful to express the distributional derivative 
of u. If Varg?„(w; [0,r]) < +00, we can introduce the distributional derivative jJ-oij^^u ■ — s' of s, 
which is a finite positive measure satisfying 

MKi,n([a,6]) = s(&) -s(a), [ C(i) d/X3j,,„(t) = - / Cit)s{t) dt for every C £ €[5(0, T). 

Jo Jo 

Notice that a singleton {t} has strictly positive measure if and only if t G J„; more precisely 

M:Ri,n(W) = o ifteC„ = (o,r)\J„, 
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with obvious modification for t = 0, T. As a general fact we have the representation formula (recall 
that t is the inverse of s) 

T S 

^# ("^' I (0,8)) =^3^1.- i-e. ^ CWdM3Ji,u(i) C(t(s))ds, (6.29) 

for every bounded Borel function ( : [0,T] — > R. Since t is injective in := t^^(C„) C (0,5), a 
Borel subset A of C„ is .2'^-negligiblc if and only if t{A) has yLtXi,u-nieasure 0. Therefore, as the 
derivatives t, u are Borel functions defined up to a ^^-negligible subset of (0, S), the compositions 
t o s, u o s are well defined in C„. The next lemma shows that they play an important role. 

Proposition 6.11. The Lebesgue measure and the vector measure u'^^ = u'^ + u'^ are 

absolutely continuous w.r.t. ^oii,u, o-n-d we have 

. du' 

= t o s and - — = u o s fi^-^ „-a.e. in Cu- (6.30) 

d^DJi,u d/^3Ji,„ 

Proof. The absolute continuity of both measures is easy, since Jtf^ < n^yi^^u by (|6.29p and the total 
variation is absolutely continuous w.r.t. /X3Ji,u thanks to (|6.14p . The first identity of (|6.30p 

can be proved as in [35i Lemma 4.1]. Concerning the second one, let us set for every smooth 
function C e Cg°(0,T) 

uo ■.= Y.mHt+)~u{t^)). 

and let us observe that we have 



(Cot)'(s)u(s)ds= / C(t(s))u(s)d,s + J,(C). (6.31) 
c„ Jc,. 

Indeed; denoting by At — (at,bt) = t^^(t), t E Ju, the connected components of [0, S] \ Cu, and 
recalling that u(at) = u(t-), u(bf) = u{t+), we have 

- / (Cot)'(s)u(s)ds = - Acot)'(s)u(s)ds+^ r(C°t)'(s)u(s)ds 

"'Cu Jo jgj^ J at 

= [\it{s))u{s)ds~Yl t\Cot){s)u{s)ds + MO^ I C(t(s))u(s)d5 + J„(C). 
Jo tgj^ J at ic„ 

Therefore, there holds 

at) du'(t) = - [ at) u{t) dt = - aHs)) u{t{s)) t(s) d,s = - / c(t(s)) u{t{s)) t(s) ds 



(Cot)'(s)u(s)ds- C(t(s))u(s)ds + J„(C) 

at)ui<t))dfl3luu{t) + Ju{C)- 

where the fifth identity ensues from (|6.3ip and the last one from (|6.29p . Since 

c(t)d<„(t)= r at)du'{t)~jao 

Jo 

we conclude the second of (|6.30p . □ 

Corollary 6.12 (Integral expression for Varg^). Let u : [0,T] ^ X fulfil Varajf, (w; [0, T]) < +oo, 
let ^ be a positive finite measure such that ^ ^ and u'^^ <^ ^, and let us set 

Jmpgj^(u; [a, 6]) := A3j„(a; M(a), u(a+)) + Ajj^ (6; «(&_), ^(6)) 

+ ^Ji„{t;u{t^),u{t),u{t+)). (6-32) 
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Then, 

Var3j„(u;[a,6]) = j D?o ( (t, ^^W); (^W, ^ W) ) d/^W + JmPKo(^; K ^D- (6-33) 

Proof. Since the expression on the right-hand side is independent of the measure /i, it is not 
restrictive to choose [i = yUKi,u; by ()6.28p we have 

6-a + VarKo(?^;[a,fo]) = / D?i (x(s); x(s)) ds + if i((s(a), s(6)) \ C„) 

J(s{a),5(b))nC„ 

3?i (x(s(t)); x(s(t))) d^ + /i([a, 5] n J„) 

(a,b)nC„ 

((t,M(t));(t(s(t)),u(s(t)))) +Jmp3j^,(u; [a, 6]), 

(a,b)nC„ ^ ' 

and we conclude by (|6.30p . □ 

6.3. Total variation for BV solutions. We focus now on the particular case (|6.10p of Example 
16. 3i when the dissipation 3^ is associated with the vanishing viscosity contact potential ^p. 

Theorem 6.13 (Comparison between Var-So (w; [•, •]) and Varp^g (w; [•, •])). For every curve u S 
BV([0,T];X) and every interval [a,b] C [0,T] we have 

Vaipx{u; [a,b]) < Varso(w; [a,b]), (6.34) 

and equality holds in ()6.34p if and only if u satisfies the local stability condition dSiocD on (a,b). 
Furthermore, if Vai's,g{u; [a,b]) < +oo, then u satisfies dSiod ) on {a,b). 

Proof. Let us first notice that the jump contributions to the total variations Varp.g and Var^o are 
the same by (|6.12p . Inequality (|6.34p then follows by applying (|6.33p and observing that for /x-a.a. 

te[o,T] 

13o((t,.(0),(^(t),^(t))) -^(^W>^W;0,^0) >*o(^(t)) (6.35) 

(where we have used the notation w{t) = — D£t(u(t))), the latter inequality ensuing from (|5.9|) . 
On the other hand, in view of ()5.10|) . ()6.35p is an identity if and only if w{t) G K* for /i-a.a 
t e (0,T), i.e. if the local stabihty property jSiod holds. 

Finally, since < n, ^^{t) > for if^-a.a t £ (0,T). Therefore, on account of ([STT]) we 
conclude the last part of the statement. □ 

Corollary 6.14. A curve u : [0,T] X is a BV solution if and only if it satisfies one of the 
following (equivalent) two conditions: 

Varso(u; [to, ^i]) + £ti (^(ti)) = £t„(w(to))+ / dtE.t{u{t)) At for every 0<to<ti<T, (6.36) 

Vars„(u; [0,T]) + £t(u(T)) < £o(^i(0))+ / dt£s{u{s)) ds. (6.37) 

Jo 

Lemma 6.15. Suppose that G AC([0,T];X), e > 0, is a family pointwise converging to u as 
£ I 0, and We : ^^T] ^ X* satisfies \\we{t) +^^t{ue{t))\\x-' ^ ^ uniformly in [Q,T]. Then, 



lim inf 



(*e("s) + "^liweit))) dt > Var3o(?/; [0,r]) > Var^xi^; [0,T]). (6.38) 



Proof. Choosing a finite partition = to < < ^2 < • • • < ^jv = T of the time interval [0, T], 
Lemma 16.61 yields 

hminf / h,{u,) + ni^eit))) dt > ^ As„((t,_i,u(t,_i)), (t„ «(*,))). 
Taking the supremum of the right-hand side with respect to all partitions of [0, T], we end up with 

mm . □ 
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Wc conclude this section with the proof of Theorem 15.81 

Proof. Let (t, u) be a parametrized soKition as in the statement of the theorem. It is easy to check 
directly from definitions (|6.1ip and (|6.19|) that 

Vars„(w;[0,T])< / ^(t(s), u(s); 0, -D£t(.) (u(s))) ds 
Jo 

<£o(u(0))-£t(s)(u(S))+ / 9t£t(,)(u(s))t(s)ds 

Jo 

= £o(w(0)) - £t{u{T)) + I dt£t{uit)) dt , 

Jo 

where the second inequality ensues from (|5.22p . Thus, (|6.37l) holds, so that w is a BV solution by 
Corollary 16. 141 The converse implication follows from Proposition 16. 101 □ 
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